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Abstract 

According to the Kosterlitz-Thouless- Theory two-dimensional solid films melt 
by the unbinding of dislocation pairs. A model including quenched random 
impurities was already studied by Nelson [Phys. Rev. B 27 (1983) 2902] , 
who predicted a reentrance into the disordered phase at low temperatures and 
weak disorder. New investigations of the physically related XY-model [e.g. T. 
Nattermann et al., J. Phys. (France) 5 (1995), 565] and a work of Cha and 
Fertig [Phys. Rev. Lett. 74 (1995) 4867] refuse this reentrant melting. In this 
work we map the system onto a two-dimensional vector Coulomb gas and via 
a renormalization we derive flow equations both for the square and for the 
triangular lattice. An analysis of these flow equations shows a new behaviour 
in the low- temperature range, where the reentrance into the non-crystalline 
phase with short-range order is not found, but the crystalline phase with quasi- 
long-range order is preserved below a critical disorder strength of <f c = j^. 
Finally we estimate the influence of commensurate substrates and obtain phase 
diagrams, which show that the melting by dislocation unbinding can only be 
expected, if the lattice constant of the crystalline film is a multiple of the lattice 
constant of the substrate potential. 
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Kapitel 1 

Einleitung 



In der Festkorperphysik versteht man unter einem dreidimensionalen Kristall 
eine (fast) unendliche Wiederholung identischer Struktureinheiten, die Translationssymmetrie 
aufweist. Der Anteil der Struktureinheiten, die die Oberflache bilden soil verschwindend 
gering sein. In drei Dimensionen bilden alle Stoffe - sieht man von Helium und 
Quasikristallcn wic Al 86 Mn 14 ab - solche kristallinen Strukturen, die aber von 
eingefrorener Unordnung gestort werden konnen. 

Daher gehdren auch dreidimensionale Kristallstrukturen wie Metalle zum alltaglichen 
Erfahrungsbereich und das Schmelzen der Kristalle zu Flussigkeiten bei Temperaturerhohung 
kann unter geeignetem Druck beobachtet werden. Hierbei handelt es sich um 
einen Phascniibergang erster Ordnung, bei dem ein Sprung in der Entropie 
auftritt. Ausgelost wird er durch die energetische Anregung des Gitters. 

Obwohl sie auf Oberflachen haufig auftreten sind zweidimensionale kristalline 
Filme hingegen nicht so leicht zu beobachten, dennoch kann man das Schmelzen 
auch hier an einigen Systemen untersuchen. Dies sind unter anderen: 

• Elektronen, die in einer einzigen Schicht etwa 100 A iiber der Oberflache 
fliissigen Heliums aufgebracht werden. Dies geschieht mittels eines elektrischen 
Felds, welches senkrecht zur Oberflache steht [Q . Bei geniigend geringer 
Dichte verhalten sie sich klassisch. 

• Polystyrolkugeln mit einem Durchmesser von etwa 0,3 /urn bilden in Losungen 
mit geeignetem pH-Wert Kolloide. Bringt man diese zwischen zwei Glasplatten 
mit einem Abstand von etwa 1 /im, formen diese Kugeln ein zweidimensionales 
Dreiecksgitter, nehmen also Kristallstruktur an (|] . 

Die Phaseniibergangs-Temperatur ist dabei systemabhangig. Sie wird von der 
Gitterkonstanten und den Lame-Kocffizicntcn bcstimmt. Dies fiihrt dazu, dafi 
die erwahnte Elektronenschicht eine Schmelztemperatur von etwa 1 K hat, 
wiihrend die Schmelztemperatur der Kolloide in der Grofienordnung von 10 2 
K liegt. 

Im Gegensatz zum dreidimensionalen Fall, gibt es bei freien zweidimensionalen 
Kristallstrukturen fur T > aufgrund von Fluktuationen in den Moden der 
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niederenergetischen phononischcn Anrcgungen keine echte langreichweitige Ordnung, 

somit also audi keine beliebig langreichweitige Translationssymmetrie. Hier existiert 

eine quasi-langreichwcitige Ordnung, bei der die Korrelationsfunktion des Ordnungsparameters 

algebraisch abfallt. Auch das Schmelzen kann in zwei Dimensionen einen anderen 

Charakter haben. Neben dem Phaseniibergang erster Ordnung existiert die 

Moglichkeit eines kontinuierlichen Kosterlitz-Thouless-Phaseniibergangs, bei dem 

sich Defektpaare auflosen und freie topologische Defekte entstehen, was zur 

Auflosung der Kristallstruktur fiihrt. Die oben genannten zweidimensionalen 

Systeme zeigen diesen kontinuierlichen Schmelzprozefi. Exakte Messungen sind 

hierbei nicht einfach, dennoch bestatigen die Ergebnisse diese im folgenden 

naher erlauterte Theorie. 

Dreidimensionale Kristalle sind gegenuber topologischen Defekten wesentlicher 
stabiler als zweidimensionale, daher tritt in drei Dimensionen kein durch topologische 
Defekte induziertes Schmelzen auf || . 

Die Arbeit gliedert sich wie folgt: Nach einer einleitenden Vorstellung des Koster- 
litz-Thouless-Phaseniibergangs mit dem Hinweis auf verwandte Systeme, wird 
ein Modell fur zweidimensionale kristalline Filme mit Unordnung motiviert 
und in eine Coulombgas-Darstellung uberfuhrt. Die mit diesem Modell von 
Nelson Q gewonnenen Ergebnisse werden diskutiert. Es folgt anschlieBend eine 
Neuuntersuchung des Modells mittels eines dielektrischen Formalismus und mit 
Hilfc der Renormierungsgruppe. Nach Presentation der Ergebnisse, die sowohl 
fur das Quadrat-, als auch fur das Dreiecksgitter keinen Wiedereintritt in die 
nicht-kristalline Phase bei tiefen Temperaturen zeigen, wird abschliefiend der 
Einfhifi eines Substrats untersucht. 



1.1 Der Kosterlitz-Thouless-Phaseniibergang 

Es ist seit langerem bekannt, daB topologische Defekte in zwei Dimensionen 
interessante Phaseniibergange verursachen konnen. Einen theoretischen Ansatz 
zur Untersuchung dieser Art von Phaseniibergangen lieferten Kosterlitz und 
Thouless zu Beginn der 70er Jahre J^, |l^| . Verschiedene zweidimensionale 
Systeme mit topologischer Unordnung, besitzen einen Phaseniibergang, der mittels 
der Kosterlitz-Thouless-Theorie erklart werden kann: 

• Das sogenannte XY-Modcll ist ein zweidimensionales Spinsystem, dessen 
klassische Spins S(x)=S(cos0(x), sin0(x)) (sic stcllen hier auch den Ordnungsparameter 
dar) topologische Defekte in Form von Vortices ausbilden konnen. Die 
Spinphase 6*(x), die die Richtung des Ordnungsparameters beschreibt, 

andert sich bei Umlauf um einen Vortexkern um 2k7r, wobei k die Windungszahl 
des Vortex genannt wird. Somit ist ein Vortex ein Spinwirbel. Experimentell 
kann ein solches zweidimensionales Spinsystem auch durch dreidimensionale 
Systeme mit schwach gekoppelten Schichten realisiert werden. Der Phaseniibergang 
entspricht hier den Ubergang von der ferromagnetischen zur paramagnetischen 
Phase. 

• Ferner ist eine Anwendung der Theorie auf suprafluide Hcliumfilme mog- 
lich, da auch hier Vortices gebildet werden konnen. Der Ordnungsparameter 
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wird allerdings nicht durch cincn zwcidimcnsionaler Vektor, sondern eine 
aquivalente komplcxc Zahl dargestellt. Das suprafluidc Helium geht mittels 
des Kosterlitz-Thouless-Phaseniibergangs in den normal fliissigen Zustand 
iiber. 

• SchlieBlich fuhren topologische Defekte auch in den im folgenden behandelten 
zweidimensionalen kristallinen Filmen zu einem Phaseniibergang von der 
kristallinen in die nicht-kristallinc Phase. Hier licgen allerdings keine Vortices, 
sondern Versetzungen vor, die sich nicht allcine durch eine skalare Windungszahl 
beschreiben lassen, da sie auch eine Richtung haben. 

Anschaulich entspricht eine Versetzung im Quadratgitter dem Endpunkt einer 
zusatzlich eingefugten Reihe von Gitterbausteinen. Dies ist in Abbildung 1.1a 
dargestellt. Der Burgers- Vektor, der entsprechend der Abbildung konstruiert 
wird, gibt Richtung und Starke der Versetzung an; sind zwei Reihen eingefugt 
ist der Burgers- Vektor doppelt so grofi, wie bei einer eingefugten Reihe. Wie in 
Abbildung 1.1b gezeigt, besteht eine Versetzung auf dem Dreicksgitter in der 
Kernregion aus einem Platz mit fiinf und einem Platz mit sieben Nachbarn. Die 
Versetzungen befinden sich auf den Gitterpunkten des dualen Gitters, welches 
beim Quadratgitter cbenfalls cin Quadratgitter, beim Dreiecksgitter ein hexagonales 
Gitter ist. Es handelt sich ausschlieBlich um sogenannte "Stufen- Versetzungen" , 
da "Schrauben- Versetzungen" nur obcrhalb von zwei Dimensionen auftreten 
konnen. 




(a) (b) 

Abbildung 1.1: Elementare Versetzungen (a) im Quadrat- und (b) im Dreiecksgitter. 
Die durch dickere Linien eingezeichneten Wege sind wegen der Versetzungen 
nicht geschlossen. Die Burgers-Vektoren sind durch Pfeile dargestellt. Die Dreiecke 
markieren die Pldtze die unregelmafiige Bonds aufweisen ■ 



Die Energie einer einzelnen Versetzung mit Burgers- Vektor b divergiert logarithmisch 
mit der Systemgrofie R, es gilt E oc b 2 In R. Daher ist es bei tiefen Temperaturen 
die Wahrschcinlichkcit, cin grofies System mit einer einzelnen Versetzung vorzufmden, 
verschwindend gering. Betrachtet man den Burgers- Vektor als Vektorladung, so 
mufi das System verschwindende Gesamtladung Ebi = aufweisen. Nur dann 
bleibt die Gesamtenergie auch fur beliebig grofie Systeme endlich, da sich die 
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obigen divergenten Beitrage in diesem Fall mit Beitragen aus der Wechselwirkung 
der Versetzungen wegheben. Der Kosterlitz-Thouless-Phaseniibergang unterscheidet 
allgemein zwischen einer Tieftemperatur- und einer Hochtemperaturphase. In 
der Tieftemperaturphase sind gepaarte Versetzungen (Versetzungsdipole) 
neben den Phononen die elementaren Anregungen des Systems. Einen solchen 
Dipol fiir das Quadrat gitter zeigt Abbildung 1.2. 
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Abbildung 1.2: Versetzungsdipol auf dem Quadratgitter mit h 1 — — t>2 = — e y 



Wie Mermin 1968 gezeigt hat jl2| , wird in zweidimensionalen Kristallen die 

langreichweitige Ordnung bei T > durch Fluktuationen zerstort, die durch 

langwellige Phononen hervorgerufen werden. Es herrscht hier quasi-langreichweitige 

Ordnung vor, was heifit, daB die Translations-Korrelationsfunktion des Translationsordnungsparameters 

pa (r) zwischen den Platzen algebraisch abfallt |n|. Der Translationsordnungsparameter 

ist wie folgt definiert: 

PG (r)= e lGu W, (1.1) 

hierbei ist G ein Vektor des reziproken Gitters und u(r) ein Vektorfeld, dafi 
die Verschiebung der Kristallbausteine von ihrem Gleichgewichtsplatz mifit. Die 
algebraisch abfallende Korrelationsfunktion nimmt folgende Form an: 

G(T) = (p G (r)*p G (0))ocr-*°. (1.2) 

(...) steht fiir eine thermische Mittelung. Der Exponent r\Q ist abhangig von den 
Lame-Koefhzienten A, fi, die Konstanten des Gitters sind 

|G| 2 r(3/i + A) 
VG = 4^2/i + A) ' (L3) 

Separiert man die phononischen Systemanregungen von den Anregungen durch 
Versetzungen und transformiert den Anteil letzterer auf ein Coulombgas (Details 
werden in Kapitel 2 erlautert), so erhalt man fiir die Energie eines Dipols 

j-n T (, , , | | (bi(ri-r 2 ))(b 2 (ri-r 2 ))\ 

EocJ bib 2 ln ri -r 2 - , ,„ , (1.4) 

V |ri-r 2 | 2 J 
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J ist hierbei eine Kopplungskonstante, die von den Lame-Kocffizienten des 
Gitters abhangt: J = J(A, fi). 

Der Betrag der Burgers- Vektoren geht quadratisch in die Energie ein, es kann 
also in guter Naherung angenommen werden, daB nur Dipole mit Vektorbetragen 
|b| = a gebildet werden, da Dipole mit grofieren Vektorbetragen demnach 
wesentlich unwahrscheinlicher sind. Bei niedrigen Temperaturen ist nur eine 
geringe Anzahl an Versetzungsdipolen vorhanden, sie nimmt aber mit steigender 
Temperatur zu, wie auch die Dipolgrofie wachst. Es kommt zu Abschirmungseffekten, 
denn fur Dipole mit weit auseinander liegenden Burgers- Vektoren wird die 
Kopplung durch dazwischen liegende kleinere Dipole reduziert. Dennoch bleibt 
in der Tieftemperaturphase die Kopplung auch im Limes r — > oo endlich (r: 
Dipolgrofie) . 

In der Hochtemperaturphase hingegen werden die Abschirmungseffekte so 

stark, dafi man fur r — > oo verschwindende Kopplung hat ( J — > 0). Die Korrelationsfunktion 

fallt hier exponentiell ab (TTJ : 

G(r) oce-f. (1.5) 

£ ist die Korrelationslange. Oberhalb der Ubergangstemperatur T m = J/16ir 
losen sich die Versetzungspaare auf, in der Hochtemperaturphase liegen damit 
freie Versetzungen vor. 

Der wesentliche Vorgang beim Kosterlitz-Thouless-Phasenubergang ist also das 
Dissoziicrcn der Vcrsctzungsdipole. Es liegt kcine kristalline Struktur mehr vor, 
da freie Versetzungen dazu fuhren, daB die Ordnung in der Hochtemperaturphase 
nicht mehr quasi-langreichweitig, sondern kurzreichweitig ist. Die Korrelationslange 
£ bei diesem Phaseniibergang nimmt folgende Werte an: 



£ oc exp 



T-T, 

£ = oo T < T, 



T > T m 

(1.6) 



Hierbei ist b eine materialabhangige Konstante, der Exponent v nimmt auf dem 
Quadratgitter den Wert 0.5 und auf dem Dreiecksgitter den genaherten Wert 
0.3696 an Q. 

Wie Kosterlitz in gezeigt hat, gilt fur den singularen Anteil der freien 
Energie 



sing 



cxr 2 . (i.7) 



Mit (1.6) folgt daraus, daB sowohl die freie Energie als auch alle ihre Ableitungen 
bei T — > T m verschwindende singulare Anteile haben und daher kontinuicrlich 
sind. Somit ist der Kosterlitz-Thoulcss-Phascnubcrgangcin Ubergang unendlicher 
Ordnung. Man kann die Phasenubergangstemperatur nur bestimmen, wenn man 
die Kopplung auf grofien Langenskalen betrachtet, die bei T m verschwindet. 

Eine besondere Eigenschaft des Phaseniibergangs bei kristallinen Filmen ist, 
dafi in der Hochtemperaturphase zwar kurzreichweitige Ordnung bezuglich der 
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Translations-Korrelationsfunktion vorliegt, aber eine Orientierungs-Korrelations- 
funktion mit einem Ordnungsparameter 

tb(r) — e 4 / 9 ^ fur Quadratgitter 

(1.8) 

ip(r) — e 6 / 9 ( r ' fiir Dreiecksgitter 

hier noch quasi-langreichweitige Ordnung aufweisen kann. 0(r)mifit die Bond- 
Orientierung relativ zu einer festen Referenzachse und ist eine Funktion, die 
alleine vom Verschiebungsfeld u(r) abhangt. Durch einen zweiten Kosterlitz- 
Thouless-Ubergang bei Tj > T m erreicht man dann die isotrop fliissige Phase, in 
der sowohl beziiglich Translation als auch beziiglich Orientierung nur kurzreichweitige 
Ordnung vorliegt. Die Zwischenphase mit quasi-langreichweitiger Ordnung in 
der Orientierung, aber kurzreichweitige Ordnung beziiglich der Translation wird 
im Fall des Quadratgitters "tetratische" , im Fall des Dreiecksgitters "hexatische" 
Fliissigkristall-Phase genannt (T^| . Die hier auftretenden topologischen Defekte 
sind sogenannte Disklinationcn: Auf einem geschlossenen Weg um eine Disklination 
in der hexatischen Phase andert sich die Bond-Orientierung um ein Vielfaches 
von 60 Grad (in der tetratischen Phase 90 Grad). Unterhalb von liegen nur 
Disklinationspaare vor, die bei TJ dissoziieren. Dabei entspricht ein Disklinationsdipol 
einer freien Versetzung, weswegen man auch Versetzungsdipole als eine Verbindung 
von vier Disklinationen betrachten kann. Abbildung 1.3 verdeutlicht die Phasenabfolge. 



test 




(b) T 

Abbildung 1.3: Phasendiagramm (a) des zweistufigen Ubergangs von einem 
zweidimensionalen Kristall in eine isotrope Fliissigkeit, (b) zeigt ein spekulatives p-T- 
Phasendiagramm, in dem feste, fliissige, gasfdrmige und Fliissigkristall-Phase auftreten 
(nach Jl^|/ ). In (b) sind diinne Lmien Kosterlitz-Thouless-Ubergange, massive Linien 
Phaseniibergange erster Ordnung. 



Der in Abbildung 1.3 eingezeichnete, rein spckulative Phaseniibergang erster 
Ordnung zwischcn fester und isotrop fliissiger Phase ist im Rahmen der Kosterlitz- 
Thouless-Theorie dann moglich, wenn die Auflosung der Disklinationspaare rechnerisch 
vor der Auflosung der Versetzungspaare erfolgen soli. Dann fallen beide Ubergange 
zusammen zu einem Phaseniibergang erster Ordnung. 
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Die folgende Tabelle veranschaulicht noch einmal die Phasenabfolge und das 
Verhalten der wesentlichen Grofien in den verschicdcncn Phasen. 





kristallin 


tetratisch/hexatisch 


isotrop fltissig 


Versetzungen 


gepaart 


frei 


frei 


Disklinationen 


als Quartette 


gepaart 


frei 


Translations- 
Korrelation 


quasi-lang- 
reichweitig 


kurzreichweitig 


kurzreichweitig 


Orientierungs- 
Korrelation 


langreichweitig 


quasi-lang- 
reichweitig 


kurzreichweitig 



Abbildung 1.4 zeigt beispielhaft das beschriebene Schmelzen von der festen iiber 
die hexatische in die isotrop flussige Phase am Beispiel eines Dreiecksgitters. 

A separate downloading of this figure is possible. 

Abbildung 1.4: Schmelzen eines Dreiecksgitters von der festen Phase (A) iiber 
die hexatische Phase (B, C) in die isotrop flussige Phase (D). Schwarz eingefarbte 
Bereiche markieren Versetzungen und Disklinationen (aus Jjjj/ ). 

In Abbildung 1.4(A) befindet sich das System noch im festen Zustand. Man 
findet hier vier Versetungspaare und vier ungebundene Versetzungen, die aber 
nahe am Rand der Darstellung liegen, deren Partner sich also aufierhalb des 
Darstellungsbereichs befinden. In (B) und (C) ist die Translations-Korrelation 
schon nicht mehr quasi-langreichweitig, wahrend dieses Verhalten fur die Orientierungs- 
Korrelation noch vorzufinden ist. Hier befindet sich das System in der hexatischen 
Phase, in der freie Versetzungen und gebundene Disklinationen vorliegen. In der 
isotrop fliissigen Phase (D) existiert keinerlei Ordnung mehr. 

Im folgenden soil aber nur auf den ersten Phaseniibergang, also auf das Schmelzen 
der festen Phase eingegangen werden. 

Die Frage, wann die Kosterlitz-Thouless-Theorie anwendbar ist und wann ein 
Phaseniibergang erster Ordnung auftritt, der nicht auf einem Defektschmelzen 
beruht, ist noch nicht vollig geklart Jl(|. Allerdings ist anzunehmen, dafi der im 
weiteren untersuchte Kosterlitz-Thouless-Ubergang nur dann auftritt, wenn die 
hicrfiir notwendigen Bedingungen (beispiclswcisc cine klcinc Vcrsetzungskonzcntration) 
erfiillt sind. Starke Anharmonizitaten wie eine starke Wechselwirkung mit einem 
Substrat konnen eine weitere Quelle fur einen Schmelziibergang erster Ordnung 
sein. 

Unter realen Bedingungen muB ein kristalliner Film immer auf ein Substrat 
aufgebracht werden. Gleichung (F4) gilt in obiger Form nur fur glatte Substrate. 
In Kapitel 5 wird versucht werden, den Substrateinflufi abzuschatzen. 
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1.2 Der Einflufi von Unordnung 

Dem kristallincn Film kann zusatzlich Unordnung in Form von eingefrorenen, im 
System nicht beweglichen, Verunreinigungen hinzugefugt werden. Realisierbar 
ist dies mit Fremdpartikeln, die Kristallplatzc besetzen und so Verzerrungen 
der Kristallstruktur hervorrufen. Da die Verteilung der Verunreinigungen in 
realen Systemen als zufallig angenommen werden kann, ist auch die lokale 
Konzentration der Verunreinigungen zufallig um einen Mittelwert verteilt. Mifit 
die Funktion Sc(r) die Abweichungen von diesem Mittelwert, so kann man 
annehmen, dafi eine gegebene Konfiguration Sc(r) mit der Wahrscheinlichkeit 

V [5c{v)] oc exp f-i- J d 2 r [<5c(r)] 2 ^ (1.9) 

auftritt B Somit sind die lokalen Abweichungen i5c(r) unkorreliert und gaufiverteilt 
mit der Varianz a. 

[Sc(t)Sc{t')] d = aS(r - r'). (1.10) 

Hierbei bedeutet [...]d eine Mittelung iiber alle moglichen Unordnungskonfigurationen. 
Die Starke der Unordnung ist gleichzusetzen mit der Starke der Fluktuationcn 
und wird durch die Varianz a ausgedriickt, die einer eingefrorenen Temperatur 
entspricht. 

Das Schmelzen von solchen kristallinen Filmen mit eingefrorener Unordnung 
wurde bereits 1983 von Nelson untersucht [gj . Abbildung 1.5 gibt das von ihm 
gefundene Phasendiagramm wieder. Auffallig ist, dafi die quasi-langreichweitige 
Ordnung, die nur unterhalb einer kritischen Unordnungsstarke a c existiert, fur 
niedrige Temperaturen T < T_(cr) wieder zerstort wird. Es ffndet demnach ein 
Wiedereintritt in die ungeordnete Fliissigkristall-Phase (hexatische, tetratische 
Phase) statt. 




Abbildung 1.4: Phasendiagramm fur kristalline Filme mit Unordnung nach Nelson 
1983. 



Ein analoges Phasendiagramm haben Rubinstein, Shraiman und Nelson 1983 
auch fur das physikalisch verwandte XY-Modell mit Phasenfrustration gefunden 



17 1. Simulationen von Josephson-Gittern zeigten allerdings keinen Wiedereintritt 



in die ungeordnete Phase |18j, 191. Ozeki und Nishimori konnten schliefilich zeigen, 
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daB ein soldier Wiedereintritt bei den vorliegenden Systemen unmoglich ist p0[ 

Neuere Arbeiten von Nattermann, Scheidl, Korshunov, Li und Tang ||, [|, 
|[ [| finden dann auch einen anderen Verlauf der Phasengrenze fur das XY- 
Modell mit Phasenfrustration. Sie argumentieren, dafi der in JT^ ] gefundene 
Wiedereintritt in die ungeordnete Phase nur ein Artefakt der Naherung fur 
klcine Vortexfugazitaten ist, die im Temperaturbereich T < T* falsch wird, da 
aufgrund des Unordnungspotentails hier hohe Vortexdichten auftreten konncn 
und somit die Wechsclwirkung zwischen den Vortexdipolen wichtig wird. Eine 
Neubetrachtung liefcrt im Temperaturbereich T < T* eine der Temperaturachse 
parallele Phasengrenze bei a c , unterhalb der die quasi-langreichweitige Ordnung 
auch bei beliebig kleinen Temperaturen erhalten bleibt (vgl. Abbildung 1.6). 
Neuere Simulationen konnten dieses Ergebnis bestatigen Q . 




Abbildung 1.5: Modifiziertes Phasendiagramm fur 
Phasenfrustration. 



XY-Modell mit 



Aufbauend auf den in || |], [|| entwickelten Methoden sollen in dieser Arbeit 
modifizierte Flufigleichungcn und Phasendiagramme auch fur kristallinc Filme 
mit eingefrorener Unordnung gefunden werden. 
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Kapitel 2 

Das Modell 



2.1 Herleitung des Modells 

In diesem Abschnitt wird ein Modell fur zweidimcnsionale kristallinc Filme mit 
eingefrorenen Verunrcinigungcn hcrgeleitet. Mit Hilfe der Elasitizitatstheorie, 
die die Mechanik fester Korper beschreibt, erhalt man eine Hamiltonfunktion 
fur das System. Hierbei werden die Verunreinigungen als lokale Quellen von 
Verzerrungen behandelt. Solche Verunreinigungen konnen beispielsweise Unebenheiten 
des Substrats, oder aber Fremdelemente im Gitter des kristallinen Films sein. Es 
handelt sich also um positionelle oder kompositionelle Unordnung. Abbildung 
2.1 zeigt ein einzelnes grofieres Atom in einem Gitter gebildet aus kleineren 
Partikcln. Wie man unschwer erkennt, resulticrcn Dchnungcn und Kompressionen 
des Gitters. 

Eine solche Verunreinigung soil eingefroren sein, was heifit, dafi sie zwar durch 
langwellige Phononen angeregt werden kann, aber dennoch nicht in der Lage 
ist, mit anderen Gitterbausteinen den Platz zu wechseln. Dies erfordet kleine 
Verunreinigungs-Diffusionskonstantenbei tiefen Temperaturen. Die Versetzungcn 
hingegen sollen sich im thermischen Gleichgewicht befinden, welches sie mittels 
Gleiten durch den Kristall erreichen. Wir betrachten also Messungen nach Zeiten, 
die grofier sind als die Relaxationszeiten der Versetzungen, die aber wesentlich 
kleiner sind als die Relaxationszeiten der Verunrcinigungcn, dcrcn Frcihcitsgradc 
eingefroren sind und deren Konfiguration wahrend der Messung somit unveraudert 
bleibt. 



2.1.1 Kristalliner Film ohne Unordnung 

Der zweidimensionale Kristall wird unter der Einwirkung innerer und aufierer 
Krafte bis zu einem gewissen Grade deformiert. So erhohen auch Versetzungen 
die elastische Energie des Kristalls. In der ublichen Kontinuumsnaherung |^2|, 
p3{ wird die Deformation durch den Verschiebungsvektor u dargestellt. Ist r 
die Position eines Kristallpunktes ohne Deformation und r' die Position mit 
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Deformation, so gilt 

u = r'-r. (2.1) 

u(r) ist damit das Verschiebungsfeld. Betrachtet man den Abstand zweier beliebiger, 
infinitesimal voneinander entfernten Punkte vor der Deformation dl und nach 
dcr Deformation dl' dann hat man mit Summcnkonvcntion 



dl' 2 = dx' 2 = {d Xl + dm) 2 
— dl 2 + 2uikdxidxk 



(2.2) 



dabei ist 



/ r N = 1 f dujijc) | du k {r) | duijr) dui(r) \ 
2 V dxk dxi dxi dxk J 

Der Tensor Uik ist der symmetrische Verzerrungstensor. Die auftretenden Deformationen 
sind stets klein in dem Sinnc, daB die Anderung eines Abstandes im Kristall 
immer kleiner ist als der Abstand selbst. Daraus resultierend sind auch die 
Verschiebungen Uj und deren Ableitungen immer klein, so daB es geniigt im 
Verzerrungstensor nur Ableitungen erster Ordnung zu beriicksichtigen. 

1 / du t (r) du k (r) \ 

^ (r) = 2 b^r + ^r • (2 - 4) 



Da eine konstante Verschiebung keine Anderung der Energie im Kristall hcrvorruft 
und da der Verzerrungtensor klein ist, kann man die Hamiltonfunktion nach dem 
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Verzerrungstensor entwickeln. Asymmetrische Beitrage niedriger Ordnung kann 
die Hamiltonfunktion nicht enthalten, wenn sie rotationsinvariant bleiben soil. 
Ferner kann die Entwicklung keine in Uik linearen Glieder enthalten, soil doch 
fur u.ik = ein Minimum vorliegen. Somit ist die allgemeinste Form, die die 
Hamiltonfunktion in niedrigster Ordnung annehmen kann: 

H = H + J d 2 rK % j k iUi u k i (2.5) 

Betrachtet man nur den von der Deformation abhangigen elastischen Beitrag 
in H und beriicksichtigt die in einem kontinuierlichen System vorliegenden 
Invarianzen, so bleibt 

H el = lJ d 2 r (2/m? fe (r) + Au£(r)) . (2.6) 

mit den beiden unabhangigen Lame-Koeffizienten A und fi. Wie Landau und 
Lifshitz 1 23 gezeigt haben, gilt diese Kontinuumsform des Hamiltonians exakt 
auch fur das zweidimensionale Dreiecksgitter. Fur das Quadratgitter findet man 
drei unabhangige Koeffizienten, die obige Form der Hamiltonfunktion kann aber 
als Naherung benutzt werden. 

Da die inneren Krafte einer Kristallflache durch von auBen, iiber den Flachenrand 
wirkende Krafte kompensiert werden, muB fiir jede Komponente des Kraftvektors 
gelten 



Fi 



rf 2 ^(r) - / dx a l3 {r). (2.7) 
Jan 



Nach den Integrationsregeln ist hierbei /; die Divergenz des symmetrischen 
Spannungstensors atj 

fi = (2.8) 

OXk 

Betrachtet man nun eine kleine Anderung des Verschiebungfelds u(r) um den 
Wert <5u(r), so ist die hieraus resultierende Arbeit: 



SW= tfrfiSui (2.9) 
Partielle Integration und Ausnutzung der Symmetrie des Tensors liefert 

flT-i/A^^U^) (2.0, 

oder 

6W = -a ik (r)6u ik (r). (2.11) 



Der Spannungstensor ist also gegeben durch 

dE 
du 



^=(1^) , (2.12) 



ik J s 
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was ein weiteres Argument fur die Nicht-Existenz von linearen Gliedern des 
Verzerrungstensors in dcr Hamiltonfunktion licfcrt, konnten diese doch bei Ui k — 
zu einer Spannung 7^ fiihren. Man erhalt nun fiir die Hamiltonfunktion 
den folgenden Ausdruck || : 

Hd = \\ d 2 ra lk {v)u tk {r), (2.13) 



2 . 

mit (wegen fl2.6Q) 

<j ik (r) = 2fiu ik (r) + XSikUjjir). (2.14) 

Diese Hamiltonfunktion, die das System ohne den EinfiuB von Unordnung beschreibt, 
wurde bereits in mehreren Arbeiten ausfuhrlich diskutiert O, |l4|, Q . 



2.1.2 Beschreibung der Unordnung 

Wie Abbildung 2.1 zeigt, sind die im Rahmen dieses Modells betrachteten 
Verunreinigungen fest lokalisierte Quellcn von Dchnungen oder Kompressionen 
des Gitters, die spharische Symmetrie beziiglich Verunreinigungspositionen ri 
aufweisen. Da das durch eine Verunreinigung erzeugte Verschiebungsfeld u lmp (r) 
ferner der Divergenzbedingung Vu lmp (r) oc <5(r — ri) geniigen mufi, hat es in 
Analogie zur Elektrodynamik in zwei Dimensionen die Form 

u imp (r) = 7^ = 7 Vln|r|. (2.15) 

Dabei mifit die Konstante 7 die "Starke" des Defekts. Eine solche Verunreinigung 
bewirkt eine Anderung Qq der Kristallfilmnache, die der Defektstarke proportional 
sein mufi und eine positiven oder negativen Wert annehmen kann. Nach p5f ist 

il = 2^7^^. (2.16) 
/i + A 



Die hieraus resultierenden Korperkrafte lassen sich mit Hilfe von und 
berechnen. 

f (r) = uV 2 u im P(r) + (A + u)V(Vu im P(r)) 

(2.17) 

= n VS(r-Ti)(fi + X). 
Die Wechselwirkungsenergie einer Verunreinigung mit einem externen Verschiebungsfeld 



u(r) erhalt man hieraus mittels (2J)). Demnach nimmt die Hamiltonfunktion 
fiir die Wechselwirkung zwischen der Unordnung und dem Verschiebungsfeld 
des Systems die folgende Form an J2|, Q 

H D = -fi (fi + A) J d 2 rc{r)u kk {r), (2.18) 

wobei c(r) = '^2 i S(r — ri) iiber alle Verunreinigungspositionen summiert. Auf 
geniigend grofien Wellenlangen kann man die Diskretheit vernachlassigen und 
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c(r) als glatte Funktion behandeln. Die lokale Konzentration der Verunreinigungen 
ist nun zufallig um einen Mittelwert Co = [c(r)],D verteilt. Es gilt 

c(r) = c + Sc(v). (2.19) 



Das Verhalten von 5c(r) ist bereits in (1.9) und ( 1.10 ) beschrieben worden. Der 



zu Co proportionale Term von Hp verhalt sich wie ein hydrostatischer Druck, 
der durch eine Dehnung des Gitters kompensiert werden kann Q . Somit ist nur 
der von Sc(r) abhangige Anteil des Hamiltonians von Intcresse. 

Die Gesamthamiltonfunktion des Systems mit Unordnung lautet demnach 

H = H e i + He, 

= \ J <fr[2^ul{v) + A^-(r)] - f7 ( M + A) J <fr5c(r) Ujj (r). ^ 

Ziel des nachsten Abschnitts ist es, diese Hamiltonfunktion auf eine Coulomb- 
gas-Beschreibung abzubilden. 



2.2 Abbildung auf ein Coulombgas 

Das Verschiebungsfeld u(r) kann in zwei Anteile zerlegt werden 

u(r) = v(r) + u*(r). (2.21) 

v(r) stellt den Verschiebungsanteil, der durch Phononen hervorgerufen wird 
dar. Dieser Anteil bildet keine topologischen Defekte aus. Fur alle geschlossenen 
Wege 7 im zweidimensionalen Film gilt 



j> dv(r) = 0. (2.22) 



Aus diesem Grund variiert der phononische Verzerrungstensor Vik (r) analytisch, 
ohne Singularitaten aufzuweisen. 

Die Verschiebung aufgrund von Versetzungen u* (r) hingegen erfiillt fur geschlossene 
Wege 7, die eine Versetzung mit Burgers- Vektor b einschliefien, entsprechend 
Abbildung 1.1 die Bedingung 



^"du*(r)=b. (2.23) 



Ferner wahlen wir u* (r) so, dafi ein System fur welches u(r) = u* (r) gilt, sich bei 
gegebener Versetzungskonfiguration in einem lokalen Energieminimum bcfindct. 
Somit diirfen in einem solchen System keine inneren Krafte vorliegen. Wegen 



(2.8) muB die Gleichgewichtsbedingung 

Oct 



ik 



dxu 



= 0, (2.24) 
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iiberall, auBer in den Versetzungskernen, erfiillt sein. Hierbei ist 
o-*fc( r ) = 2 A"4(r) + XSikUjjir). 



(2.25) 



Setzt man diese Zerlegung in die Hamiltonfunktion ein, so erhalt man einen 
phononischen und einen Versetzungsanteil, sowie einen Term der die Wechselwirkung 
zwischen Phononen und Versetzungen beschreibt 



wobei 



H e l — H p hon + Hver + HwW j (2.26) 

H phon = i/d 2 r(2 M ^ fe (r) + A^(r)) 

H Ver = i/d 2 r(2 M < fe 2 (r)+A< t 2 (r)) (2.27) 

Hww = ^Jd 2 r (2fjLv ik (r)u* ik (r) + Xvu(r)u^(r)) . 

Der Wechselwirkungsantcil nimmt nach explizitem Einsetzen des Tensors Vn- 
folgende Form an 



1 f , 2 / / dvi dvk\ dvi 
Hww = - J d r f 2 M f — + ^; j u ifc + 2A— Uji 



(2.28) 



Dieser Ausdruck kann partiell integriert werden, wobei keine Randterme auftrcten, 
da wir annehmen, da8 das Hinzufiigen der Versetzungen die Form des Systemrandes 
nicht verandcrt. Daher 



WW 



d 2 r 



dxi 



— — + v k — — 

dxi dxk 



(2.29) 



ox k 



Vk- 



dxi 



Aufgrund von ( 2.24 ) ergibt sich Hww = 0, die Hamiltonfunktion zerfallt also 
in einen rein phononischen Anteil und einen Anteil, der nur die Wechselwirkung 
der Versetzungen untereinander beschreibt. Da wir im folgenden am Verhaltcn 
der Versetzungen interessiert sind, bctrachtcn wir nur Hver- 



2.2.1 Abbildung des Systems ohne Unordnung 



Fur ein System ohne Unordnung nimmt die zu untersuchende Hamiltonfunktion 
nach ( 2.13 ) folgende Form an 

1 



H = 



d ra* k (v)u* k {v). 



(2.30) 



Wegen ( 2.24 ) gelten die folgenden beiden Gleichungen iiberall, aufier in den 
Versetzungskernen 

dot 



dx 



d(T U , da *yy 
dx 



dy 



dy 







= 



(2.31) 
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Unter Beriicksichtigung der Symmetric von a* k , ist die folgende Darstellung von 
a* k die allgemeine Form der Losung dieser beiden Gleichungen pj| : 

o-S b (r)=eye w V J Vix(r) (2.32) 

Die Funktion x nennt man Verzerrungsfunktion und hat im Potential <3? ihr 
elektrodynamisches Analogon. Aufgrund der Beziehung zwischen a* k und dem 
Verzerrungstensor u* k findet man 

<fc(r) = i-eye W ViVa(r) - A + V 2 X (r)<5 tfc . (2.33) 



Fiihrt man eine zweidimensionale vektorielle Versetzungsdichte b(r) = 
$^ a b a £(r — r a ) ein, welche uber alle Versetzungspositionen r Q summiert, so 
lafit sich ( 2.23| ) schreiben als 



j> du*(r) = j> Wji(r)dxj = / e k j\7 k Wji(r)d 2 r = (2.34) 



Hierbei ist tUji = V^u* , und es wird nur iiber die Burgers- Vektoren h a summiert, 
die von dem Weg 7 eingeschlossen werden. Driickt man diese Gleichung mittels 
der Versetzungsdichte aus, so resulticrt 



e/cjVfeWj^r) = bi(r). 



(2.35) 



Wendet man nun EjjEfeiVjV; auf beide Seiten der Gleichung (12.33 1 ) an, so ergibt 
sich mit dem Kopplungskoemzienten J = 4//(/i + A)a 2 / (2/z + A) p2|, B7j 



jV 4 x(r) = ^eijCkiVjViiWik +Wki 
= |eij£«V 'j-V 'i(w ik - w ki 
= ejiVjbiir) =: r?(r). 



(2.36) 



Der erste Term der zweiten Zeile entfallt, da wir ein System ohne freie Disklinationen 
betrachten. Die oben defmierte Funktion r/(r) nennt sich Quellenfunktion und 
beschreibt die Verteilung der Versetzungen. Sie entspricht der elektrodynamischen 
Ladungsdichte p. Der Vergleich mit der Elektrodynamik, wo V 2 $(r) — — 47rp(r) 
gilt, wahrend hier V 4 x(r) = Jri(r) ist, zeigt, dafi hier doppelt so oft abgeleitet 
wird. Dies ist notwendig, um den Vektorcharakter der betracheten Ladungen zu 
beriicksichtigen. 

Gelost wird V 4 x(r) = Jv( r ) m r |r| ^ a (die Gitterkonstante a wird im folgenden 
bis einschlieBlich Kapitel 4 immer gleich 1 gesetzt, gefordert ist also |r| ^> 1), 
mit Hilfe der Green-Funktion 



,g(r) = ^r 2 (m|r|+C) 



(2.37) 



die 



V 4 5 (r) = — V 2 (ln|r|+C + l)=£(r) 



(2.38) 
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crfiillt. Denn setzt man fur die Verzerrungsfunktion 

X {t)=J J d 2 r'r,{r')g{v-v') (2.39) 

an, so gilt ( 2.3G| ). 

(2.3C) liefert mit ( p. 32 ) und ( 2.33| ) nach zweimaliger partieller Integration folgenden 
Ausdruck fiir die elastische Hamiltonfunktion H 



(2.40) 



Die nicht beriicksichtigten Randterme liefern Beitrage die proportional zu In R 
sind, also mit der Systemgrofie divergieren. Sie entfallen allerdings, wenn das 
System Ladungsneutralitat aufweist. Um die Energie endlich zu halten, mufi 
daher 



5> Q = o. 



(2.41) 



crfiillt sein. 

Eine etwas aufwendigere Rechnung unter Verwendung der Ergebnisse fiir x und 
r\ (vgl. Anhang A.l) fiihrt schliefilich zu der gesuchten Coulombgas-Darstcllung 
der Hamiltonfunktion fiir das Problem ohne Unordnung (r aj g :— r a — rp) |l4|, 

H * = -T- E ( h ^ ln M - ^^^A +E c Y: |b Q | 2 (2.42) 



Auffallig ist hierbei der von der iiblichen Elektrodynamik abweichende zweite 
Term, der sich aus dem Vektorcharakter der hier betrachteten Ladungen ergibt. 
Die Energie ist nicht nur vom Abstand der Ladungen voneinander, sondern auch 
von deren Stellung zueinander abhangig. Dieser zweite Term nimmt beispielsweise 
fur ein Paar aus antiparallclcn Burgers- Vektoren gleichen Betrags die Form 
— -^b 2 cos 2 9 an (vgl. Abbildung 2.2). Hier ist also bei gleichem Abstand ein 
Paar, bei dem die Burgers- Vektoren parallel zum Abstandsvektor stehen energetisch 
giinstiger, als ein Paar, bei dem die Burgers- Vektoren senkrecht zum Abstandsvektor 
sind. 

Allerdings ist der Wert dieses zweiten Terms unabhangig vom Betrag des Abstandes 
der Ladungen, er wird also bei sehr groBen Abstanden klein gegenuber dem 
ersten Term. Der dritte Term ist ein Ausdruck fiir die Energie der Versetzungskerne 
(vgl. Anhang A.l) und zeigt auch keine Abstandsabhangigkeit. 



2.2.2 Abbildung des Unordnungsanteils 



Die Hamiltonfunktion fiir die Wechselwirkung zwischen einer Vcrunreinigung 



und einem auBeren Verschiebungsfeld u(r) beschreibt Glcichung (2.18). Fiir das 
externe Feld setzen wir im folgenden Formel ( p. 35 ) an. Physikalisch bedeutet 
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Abbildung 2.2: Schematische Darstellung ernes Paars aus Burgers-Vektoren 



das, dafi wir nur die Wechselwirkung zwischen der Unordnung und den Versetzungen 
beriicksichtigen, die direkte Unordnung-Unordnung- Wechselwirkung aber vernachlassigen. 
Gerechtfertigt ist dies, da diese Wechselwirkung von d or Or dnung 8c(v)8c(y') ist 
und bei Unordnungsmittelung unter Benutzung von (1.10) wegfallt B. Es gilt 
also 



1 



2A 



(/z + A)Q / d 2 rSc(r)V 2 X (r). 



(2.43) 



2n 4 ) u(A + /i) / 

Bei Verwendung der Definition der Verzerrungsfunktion \ folgt sofort 

j r i 

H D = -fi / d 2 r5c(r) ^ e zi jb ai Vj — In |r - r a |. (2.44) 



e z ij ist hier der ubliche Tensor dritter Stufe, wobei die z-Richtung, die Richtung 
senkrecht zum Kristall ist. Summiert man ( 2.44[ ) mit der elastischen Hamiltonfunkton 



( 2.42j ) , ergibt sich abschliefiend die Hamiltonfunktion, die das vollstandige Problem 



auf eine Coulombgas-Beschreibung abbildet 

i? = - 2^ b a b^ln|r a ^| 2 + E c l b « 

' J : ^ |^ C (r)^ ez[bQX(l " r " )] 



4tt u 7 (r-r a )2 



2 

(2.45) 



Offcnbar werden Versetzungen mit steigendem Betrag dcs Burgers- Vektors energetisch 
immcr ungiinstiger. Daher ist anzunehmen, dafi bei ticfen Tcmperaturen lcdiglich 
solchc Burgers-Vektoren auftreten, dcrcn Betrag der Gitterkonstantcn cntspricht, 
fur die also |b Q | = 1 gilt. Man sieht ferner, dafi sich antiparallele Burgers- 
Vektoren anziehen, wahrend parallclc Burgers-Vektoren Energie gewinnen, wcnn 
sie sich voneinander entfernen. Bei senkrecht zueinander stehenden Vektoren ist 
die Wechselwirkungsenergie unabhangig vom Betrag ihres Abstands, fallt doch 
der erste Term der Hamiltonfunktion in diesem Fall weg. Somit bilden sich bei 
geniigend tiefen Temperaturen und nicht zu starker Unordnung Dipolpaare, freie 
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Versetzungen werden erst bei hoheren Temperaturen oder starker Unordung 
energetisch moglich sein. 

Da bei tiefen Temperaturen und geringer Unordnung auch Ladungsneutralitat 
(S^a = 0) vorliegt, ist es moglich allc Versetzungen eindeutig zu Dipolen 
zu paaren. Dabei sind Dipole stets Versetzungspaare, deren Burgers- Vektoren 
antiparallcl zucinander stehen. Wegen |b Q | = 1 sind Dipole somit auch ladungsneutral. 
Zur Bildung des ersten Dipols betrachtet man jede Versetzung und sucht die 
dazu nachstliegende Versetzung mit antiparallelem Burgers- Vektor. Jenes der 
moglichen Paare mit der geringsten Separation bildet den ersten Dipol. Zur 
Bildung des zweiten Paares geht man unter Vernachlassigung der beiden Versetzungen, 
die den ersten Dipol bilden, genauso vor. Nach genugend vielen Schritten hat 
man so alle Versetzungen eindeutig gepaart. 

Die Wechselwirkung zwischen zwei Versetzungen eines Dipols mit Separation 
r wird durch Dipole mit kleiner Separation abgeschirmt. Dieser Effekt kann 
beriicksichtigt wcrdcn, in dcm man allc Dipole mit einer Separation kleiner als 
r nicht mehr als Ladungen, sondern nur noch als polarisierbares Medium auffafit. 
Dieses Vorgehen fuhrt zu einer Renormierung der physikalischen Grofien wie der 
Kopplungskonstanten, deren Wert dann von r abhangt. 



2.3 Bisherige Resultate 

2.3.1 Wiedereintritt in die ungeordnete Phase 

Die im letzten Abschnitt hergeleitete Hamiltonhmktion fur das ungeordnete 
System ( 2.4fj ) wurde bereits 1983 von Nelson untersucht j^j . Das wesentliche 



Resultat dieser Arbeit waren die folgenden FluBgleichungen, die mittels des 
Replikatricks und Renormierung, sowie mittels einer Naherung fur kleine Werte 
der Fugazitat y hergeleitet wurden: 

dJ 3 J 2 2 JT- S f ( fJT-aJ 2 \ (JT-dJ 2 \\ 
dy ( JT - a J" 1 



+ 0{y 2 ) (2.47) 

y ovi ± ~ i 

% = (2-48) 



Es wird also eine Renormierung der Kopplungskonstante J(l) und der Fugazitat 
y(l) gefunden, wahrend die Unordnungsstarke a unrenormiert bleibt. Vom Gittertyp 



hangt der Beitrag des zweiten Summanden in (2.47) ab, der hier mit 0(y ) 



bezeichnet ist. Die sogenannte Fugazitat y besitzt folgenden Ausgangswert 



yo = y(r = 1) = e t = e . (2.49) 

Auf grofieren Langenskalen r := e l ist die renormierte Fugazitat y = e~ Ec ^ r ^ T , 
wobei die renormierte Energie Ec(r) eines Versetzungskerns seiner freien Energie 
Fq = Ec' — TSc auf einer Flache der Grofie r 2 entspricht. Somit nimmt y 2 Werte 
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an, die proportional sind zur Wahrscheinlichkeit Versetzungsdipolc zu findcn 
und proportional zu r A sind, der Zahl der durch Rcnormierung verlorengegangenen 
moglichen Dipolpostionen (vgl. hierzu Kapitel 4) . Ferner wurde eine neudefinierte 
Unordnungsstarke a — <jQq und die ublichen modifizierten Besselfunktioncn 
Iq und 1 1 verwendet. Fiir a — > stimmen obige FluBgleichungen mit jenen 
Gleichungen iiberein, die fiir Systeme ohne Unordnung gefunden wurden 13, I3j 



Eine numerische Integration der FluBgleichungen licfert die in Abbildung 2.3 
gezeigte Darstcllung des Hamiltonianflusses. 




Die Anfangsbedingungen sind durch die i/q darstellende gestrichelte Linie gekennzeichnet. 
Auf der Achse y — befinden sich zwei ausgezeichnete Punkte KZ 1 und KZ 1 



± 32tt 



1 ± (1 - 64ttct) 1/2 . (2.50) 



Nur fiir KZ 1 < R- 1 < KZ 1 geht der Flufi zur y = 0-Linic hin, in alien anderen 
Fallen fliefit er von dieser Linie weg. Die spezielle, dicker eingezcichnctc, Linie, 
die die Fixlinie y = bei Kq 1 verlafit und bei KZ 1 wieder in die Fixlinie 
einmiindet, separiert zwei Phasen. Man nennt sie daher auch Separatrix. Im 
schraffierten Bereich liegt die feste Phase, da hier der FluB immer bei y = und 
endlichen Werten von K einmiindet. Wegen der oben erwahnten Proportionalitat 
zwischen der Fugazitat y und der Dipolwahrscheinlichkeitsdichte gibt es auf 
grofien Skalen offenbar kcinc Dipolc mehr, alle Versetzungen sind zu Dipolcn 
kleinerer Separation gebunden. Im nicht-schraffierten Bereich strebt der FluB zu 
grofien y- Werten; man hat also selbst fiir r — > oo eine hohe Dipolwahrscheinlichkcit. 
Ein Versetzungsdipol mit Separation r — > oo ist aber gleichbedeutend mit der 
Existenz freier Versetzungen, somit liegt hier die nicht-gcordnctc Phase vor. 

Die feste Phase existiert offenbar nur fiir Unordnungsstarken, die kleiner sind 
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als eine kritische Starke a c 

*<-a c =±-, (2.51) 

denn bei a c fallen KZ l und K+ 1 zusammen. Die realen Phaseniibergangs- 
Temperaturen ergeben sich aus den Schnittpunkten der Anfangsbedingungs- 
Linie mit der Separatrix. Im Limes Ec — > oo fallt die obere kritische Temperatur 
dann mit KT 1 zusammen. Deutlich wird, dafi auch im Falle sehr geringer Unordnung 
a > bei geniigend tiefen Temperaturen ein Wiedereintritt in die ungeordnete 
Phase stattfmdet. Das in Abbildung 1.5 dargestellte Phasendiagramm resultiert. 



2.3.2 Scheitern der Naherung fiir kleine Fugazitaten 



Nelson hat in obiger Arbeit M fiir die Wahrscheinlichkeit, einen Versetzungsdipol 
vorzufinden, durchweg eine bosonische Boltzmann - Verteilung angesetzt. Ein 
Versetzungsdipol der Separation r = \r± — mit einem Winkel 9 zwischen 
Abstandsvektor und Burgers- Vektoren (vgl. Abbildung 2.2) soil somit eine Wahrscheinlichkeit 



P{r, 9) oc 



exp 



T 



(2.52) 



D 



habcn. [...]d bedeutet wiederum eine Mittelung iiber alle Unordnungskonfigurationen, 
wahrend Ei(r,0) die Energie eines solchen Versetzungsdipols mit |b| = 1 ist. 
Nach (05} gilt 



E b (r,6) = — 6 2 (lnr-cos 2 

47T 



2b 2 E c + V{Yi) + V(r 2 ), 



(2.53) 



mit 



T7-/ \ J n f , 2 x i ^ I b a x (r- r Q )] 

= — £l ( d 2 rSc(r)(e z x b Q )Vln|r — r a | 
4vr J 



(2.54) 



Da die Energiekosten fiir Versetzungen mit oc b 2 ansteigen, wahrend ein moglicher 
Energiegewinn durch Wechselwirkung mit der Unordnung nur oc b sein kann, 
ist die Wahrscheinlichkeit Versetzungen mit einem Burgers - Vek tor |b| > 1 zu 
erzeugen verschwindend gering und wird im folgenden wie in (2.52) vernachlassigt. 



Es ist zu beachten, dafi die letzten beiden Terme der Unordnungs- Wechsel- 
wirkung in ( 2.53| ) V{y 1 )+V(v 2 ) negative Werte annehmen kSnnen und sogar den 
gesamten Ausdruck E\ (r, 9) negativ machen konncn. Bei sehr tiefen Temperaturen 
dominicrcn solchc Konfigurationcn, hier ist also in den moisten Fallen E\ (r, 9) 3> 
T nicht mehr gewahrleistet. Folglich ergeben sich hier auch sehr grofic Dipolwahrschcinlichkeiten 
P(r, 9) 3> 1, was dazu fiihrt, dafi Wechselwirkung zwischen den Versetzungspaaren 
beriicksichtigt werden miissen, man also kein verdiinntes System mehr vorliegen 
hat. Insbesondere wird die in B durchgefiihrte Entwicklung nach kleinen 
Potenzen der Fugazitat unhaltbar, ist y doch, wie bereits erwahnt, proportional 
zur Dipolwahrscheinlichkeit. 



22 



2.3. BISHERIGE RESULTATE 



Dafi dicsc Problcmatik zu cincm unphysikalischen Vcrhalten fiihrt, sollen die 



folgenden Abschatzungen zeigen. Entsprechend (2.47) gilt 



i^ 2 ^ + ^j^°^ ^ 

Die y-Proportionalitat dcs fctt gcdruckten Terms ist fur grofie Temperaturen 
plausibel, wird doch der EinfluB des Unordnungspotentials durch thermische 
Fluktuationen verringert. Im Fall a > und T — > jedoch fiihrt dieser Term 
zu einem unbegrenzt star ken Anwachsen von -# und damit, wie in Abschnitt 
2.3.1 beschrieben, zu einem Wiedereintritt in die ungeordnete Phase. Der 0(y 2 )- 
Term, der fur das quadratische Gitter null ist, nimmt im Fall des Dreiecksgitters 
lediglich kleine positive Werte an und kann daher bei dieser Betrachtung iibergangen 
werdcn. 

Versctzungsdipole mit einer Grofie zwischen r und r + dr tragen mit dFy zur 
freien Energie pro Einhcitsflache bei. Fur dr — ► kann dieser Beitrag nach 
Kosterlitz mittels 

f f 27T 

dF v fa -Trdr- J d6 [\nZ] D (2.56) 

angenahert werden (vgl. auch Scheidl || und Tang ^ ). Hierbei ist Z die 
lokalc Zustandssumme. Sie zahlt bei gegebenen Versetzungspositionen, deren 
Separation im Intervall [r, r + dr] liegt, die moglichen Dipolc. Im Falle des 
Quadratgitters sind vier Einstellmoglichkeiten der Dipole zu berucksichtigen, 
beim Dreiecksgitter gibt es sechs Moglichkeiten. 

Eine ausfuhrliche Diskussion der Einstellmoglichkeiten und orientierten Dipolwahrscheinlichkeiten 
P(r,6) erfolgt in Abschnitt 3.2 (vgl. insbesondere Abbildung 3.2). Nimmt man 
fur das Quadratgitter analog zu Nelson kleine Fugazitaten und damit kleine 
Dipolwahrscheinlichkeiten an, so hndet man wegen ln(l + x) pa x fur i<1: 

T ( 2t: it 3tt 

dF v =-~r d6rdr[P{r, 6) + P(r, 6+-) + P{r, 9 + tt) + P(r, 6 + —)} 
4 Jo 2 2 

p2tt 

= -Trdr I d6P(r,6) 
Jo 



(2.57) 

Setzt man eine Boltzmann-verteilte Wahrscheinlichkeit P{r,8) oc r _3 ^( 1_ ~) 
e s£r(i-#)<™ 2 o anj so fo i g t 



dF v oc -nTrdrr-^^-^e^i 1 - 4 ?)^ ( — ( 1 - —) 



(2.58) 



Unter Bcnutzung der von Nelson verwendeten Relation fur die Fugazitat y 2 (r) = 
ergibt sich fur das Skalenverhalten der freien Energie 
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Im Fall des Dreieckgitters findet man in niedrigster Ordnung, abgesehen von 
einem anderem konstanten Vorfaktor, dasselbe Resultat. Formal wurde dieses 
Ergebnis bereits von Kosterlitz |l0| fiir das XY-Modell ohne Unordnung und 
von Tang [|| und Scheidl || fiir den Bereich hoher Temperaturen des XY- 
Modells mit Unordnung gefunden, ihre Resultate unterscheiden sich lediglich 
durch Faktoren. 

Mittels der Beziehung = —-§f^j- ist es nun moglich das Skalenverhalten 
der Entropie zu untersuchen. Hierbei ist es nicht notig die beiden hinteren 



Faktoren aus Gleichung (2.59) zu betrachten, da sie auf grofien Langenskalen 
lnr ^> tr und lnr ^> ^f- (nur hier gelten auch die verwendeten Relationen 
fiir z.B. die Wahrscheinlichkeit) nicht zum physikalischen Verhalten beitragen. 
Somit wiirde eine Beriicksichtigung dieser Faktoren die Rechnung uniibersichtlicher 
machen, aber das Ergebnis nicht beeinflussen. 

Wir berechnen also lediglich: 

^-(Ty 2 ) =y 2 + T^-y 2 

dT f (2.60) 

=y2+ 4^ (T - 2 " J) ^ 

Hieraus folgt, dafi die Entropie im Falle T < T* = 2 J a und fiir a < a c 
unbeschrankt abnimmt. Vernachlassigt wurde bisher die bei I > vorhandene 
implizite Temperaturabhangigkeit der abgeschirmten Grofie J. Wie Scheidl || 
bereits ausgearbeitet hat, konnen diese impliziten Abhangigkeiten beriicksichtigt 
werden, indem man die Grofie -^-^(Ty 2 ) berechnet. Sind sowohl -§f(Ty 2 ) als 
auch -§fji{Ty 2 ) positiv, so ist der Entropieflufi auf jeden Fall positiv, da die T- 
Abhangigkeit von J eine Konsequenz der Abschirmung ist. Sind beide Ausdruckc 
negativ, liegt auch ein negativer Entropieflufi vor. 

Anstelle von ■^■§i(Ty 2 ) kann man ebenso 

berechnen. Fiir T > T* ist der Entropieflufi somit tatsachlich nicht negativ. 
Die unbeschrankte Abnahme der Entropie im Fall T < T* = 2Ja bestatigt 
sich allerdings. Auf grofien Langenskalen fiihrt dies zu einer negativen Entropie. 
Dieses Verhalten ist offenbar eine Konsequenz der unkorrekten Naherung fiir 
kleine Fugazitaten. Es ist daher notwendig den gesamten Bereich, in dem T < 
T*(a) gilt neu zu untersuchen. Da T*(a) die Phasengrenze gerade bei a = a c 
schneidet (vgl. Abbildung 1.6) fallt der gesamte Bereich des Wiedereintritts in 
die ungeordnete Phase darunter. 



2.4 Abschatzung der Phasengrenze 

In diesem Abschnitt soil die Phasengrenze im neu zu untersuchcnden Bereich 
T < T*(a) abgeschatzt werden. 

Da ein Phaseniibergang in die ungeordnete Phase, in der Versetzungen nicht 
mehr zu Dipolen gebunden sind, mit der Moglichkeit der Entstehung freier 
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Versetzungen durch thermische Fluktuationcn glcichzusctzen ist, untersuchcn 
wir zunachst das System bei T = und schwacher Unordnung darauf, ob das 
Unordnungspotential hier die spontanc Entstchung von Versetzungen ermoglichen 
kann. Die Encrgie einer einzelnen Versetzung in einem System der Flache Rx R 
ohne Unordnung ist, wie man aus ( p35| ) berechnen kann (vgl. ^ ) 



E el = —lnR. (2.62) 



Dies ist zu vergleichen mit dem moglichcn Encrgiegewinn Ed durch die Wechselwirkung 
V(r) zwischen der freien Versetzung und der Unordnung. Diese Wechselwirkung 
ist gauBverteilt, da auch die Funktion 5c(r) eine solche Verteilung aufweist. Die 



Varianz u>(R) der Wechselwirkung lafit sich mit Formel (2.54) berechnen (vgl. 
Anhang B) 

u>(R)= [V(r a ) 2 ] D = ^-alnR. (2.63) 
Der typische Energiegewinn durch diese Wechselwirkung ist also 



E% vp) =-j(^-lnR) ' , (2.64) 



1/2 

1 1 ■ r.» 

V167T 

und die freie Energie der einzelnen Versetzung ist damit: 

J / a \V2 

F = — In i? - 2T In i? — J ( lni? . (2.65) 

8n \l67r / 

Der zweite Term resultiert aus der Entropie, die S — ln(i? 2 ) betragt, da eine 
einzelne Versetzung R 2 Platze im System einnehmen kann. Vergleicht man bei 
T — * den ersten mit dem letzten Term, so zeigt sich, dafi bei schwacher 
Unordnung a <C 1 die Entstchung freier Versetzungen nicht zu erwarten ist. 
Mochte man aber rigoros ausschliefien, dafi das Unordnungspotential die Entstchung 
freier Versetzung induzieren kann, so mufi man nicht den typischen, sondern den 
maximal moglichen Energiegewinn j^^ 111 ^ (j er Wechselwirkung untersuchen. 
Dieser ist dann geben, wenn die Versetzung ihre energetisch giinstigste Position 
im System einnimmt. 

Hierzu folgen wir der Methode von Korshunov und Nattermann [^8| und 
vernachlassigen die Korrelationen zwischen den Werten des Unordnungspotentials 
an den verschiedenen Positionen. Wir betrachten also R 2 unabhangige Variable 
V(r) die jeweils dieselbe Gaufiverteilung 

aufweisen. Die Wahrscheinlichkeit genau den Wert ^™ ax ^> a j s Maximalwert zu 
finden ist gegeben durch die Wahrscheinlichkeit, ihn an einem Platz vorzufindcn, 
multipliziert mit der Wahrscheinlichkeit, dafi allc anderen Platze eine geringere 
Energie aufweisen und multipliziert mit der Anzahl der Platze R 2 . Die Verteilung 
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von e^ 1 "-^ igt also 



P{E ( ™ ax) 



: P( E D 



dE 



(max) 
D 



-i ir-i 



dVp(V) 



dVp(V) 



R 2 



(2.67) 



Vollig analog zu |28 lafit sich das Maximum dieser Verteilung berechnen, das 
wegen der geringen Breite dieser nicht-gaufiischen Verteilung mit dem wahrscheinlichsten 
Wert der Verteilung gleich gesetzt werden kann. Man findct 



(max) 



E 



D 



V2cj In R 2 . 



(2.68) 



Einsetzen von lo(R) licfcrt fur die gemittelte freie Energie einer Versetzung, die 
sich am energetisch giinstigsten Platz des Systems befindet 



F = — (l — Vl6ird) mi?. 

87T V / 



(2.69) 



Auch diese Gleichung zeigt, dafi fur T — ► und kleine Unordnungsstarke immcr 
noch eine im Falle R — > oo unendliche Energie aufgebracht werden mufi, um eine 
freie Versetzung zu erzeugen. Erst bei Unordnungsstarken a > kann eine 
freie Versetzung die Energie des Systems absenken und daher sind ab diesem 
Wert ungebundene Versetzungen auch zu erwarten. 

Zusatzlich zu obigen Betrachtungen kann untersucht werden, bis zu welcher 
Unordungsstarke es bei T = energetisch gunstig ist, gebundene Versetzungsdipole 
zu erzeugen. Auch die Wechselwirkung eines solchen Paares mit der Unordnung 
ist gauBverteilt. Die Varianz dieser Verteilung ist wegen ( 2.53 ) fur cinen Dipol 
mit der Separation r = |r Q — r^l gegeben durch A 2 (r, 6) — \ [V{v a ) + V(rp)) 2 ] 
Bcachtet man die Translationsinvarianz im System, so crhalt man nach einiger 
Rechnung (vgl. Anhang B) fur Separationen r 1 



A 2 (r,0)- [(V(r) + V(0)f 



2[V(0)V(0) + V(v)V(0)} 1 



Pa 

^T (lnr " 



(2.70) 



cos 



hierbei ist wiederum der Winkel zwischen den Burgers- Vektoren und dem 
Abstandsvektor. 



Da die elastische Energie eines Versetzungsdipols durch 



J 

47T 



lnr 



J 

47r' 



-COS 



2 6 + 2E r 



(2.71) 



gegeben ist, kann die Erzeugnung eines solchen Dipols bei T = nur dann 
energetisch gunstig sein, wenn die Wechselwirkung mit der Unordnung die Energiekosten 
durch E e i aufhebt. Aufgrund der Definition von Eq (vgl. Anhang A.l), kann E e i 
selbst nicht negativ werden. Entsprechend Abbildung 2.4 mufi also der Betrag 
dieser Wechselwirkung V im schraffierten Bereich liegen. 
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P(v) 




"Eel V 



Abbildung 2.4: Gaufiverteilung der Unordnungswechselwirkung eines 
Versetzungsdipols 



Daher ist die Wahrschcinlichkcit cincn Vcrsctzungsdipol mit Burgersvektor b 
bei r a und mit -b bei zu finden 

f~ Eel dV v 2 
P(r,6) = / -^e"^ 

1 A 2 _*2« 
w — — e 

/ A 2 

wobei lnr » <r angenommen wurde. Setzt man hier die Ausdriicke fur A 2 und 
E e i ein, so resultiert, da man fur grofie r die Energie der Versetzungskcrnc 
vernachlassigen kann (J lnr Ec, lnr 3> 1), 

) ~ \/ i ~ 9/ , r~^e^ cos2 e . (2.73) 
V In r — cos^ 

Die Anzahl der moglichen Platze fur einen Versetzungsdipol mit grofier Separation 
r skaliert cx r 4 . Berucksichtigt man nur die wesentlichen Abhangigkeiten in r, 
so zeigt sich, dafi die Zahl der Platze P, fur die es energetisch vorteilhaft ist, 
einen Versetzungsdipol zu erzeugen, wie folgt skaliert 

p ar 4 r -4^F. (2.74) 

Auch hier stcllt man fest, dafi im Falle a > die Zahl der Dipole mit der 
Langcnskala r ansteigt, auf beliebig groficn Skalen also auch belicbig viele Dipole 
erzeugt werden. Wie in Abschnitt 2.3.1 erklart, befindet man sich somit in der 
ungeordneten Phase. Fur a < jh- hingegen findet man das in der festen Phase 



(2.72) 
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zu erwartende Verhalten. Diesc Bcrechnung wie auch die vorherige Abschatzung 
zeigen, daB die kritische Unordnungsstarke, oberhalb der keine feste Phase mit 
quasi-langrcichweitiger Translationsordnung mehr bestehen kann, bei 

** = i (2-75) 

liegt. Dies steht im Widerspruch zu der von Nelson Q gefundenen kritischen 
Unordnungsstarke von g|^. Dieser Unterschied resultiert aus einem Rechenfehler 
in |Q , auch im folgenden werden fur T > T* die von Nelson gefundenen 
Flufigleichungen mit einer um den Faktor 4 starkeren Unordnung reproduziert. 
Der Wert a c = — ^ ist ebenso von Cha und Fertig [j?J gefunden worden. Dem 
entsprechend ergibt sich fur den Wert von T* jetzt 

T* = (2.76) 



Aufgrund vorangegangener Abschatzungen ist im gesamten Bereich T < T* 
eine Phasengrenze zu erwarten, die bei a c parallel zur T-Achse verlauft (vgl. 
Abbildung 1.6). Dieser Verlauf wiirde im Gegensatz zu nicht den Ergebnissen 
von Ozeki und Nishimori [po[ widersprechen, die einen Wiedereintritt in die 
ungeordnete Phase ausschliefien konnten und zeigten, daB die Phasengrenze 
bei tiefen T parallel zur Temperatur-Achse verlaufen mufi, falls eine Phase mit 
quasi-langreichweitiger Ordnung iiberhaupt existiert. 
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Kapitel 3 

Dielektrischer Formalismus 



In Abschnitt 2.2.1 wurde der Zusammenhang zwischen der zweidimensionalen 
Elektrodynamik und dem vorliegenden Problem dargestellt. Wie in der Elektrodynamik, 
kann man bei der Rcnormierung das Verhalten auf grofieren Langenskalen r 
mittels einer dielektrischen Konstanten e(r) beschreiben, die die Polarisierbarkeit 
von Versetzungpaaren kleinerer Separation beriicksichtigt. Da die Wahrschcinlichkeit, 
gcbundene Versetzungen zu finden, in der geordneten Phasen fur r — > oo gegen 
null geht, ist hicr cin cndlichcr Wert fur e(r) zu erwarten, wahrend in der 
ungeordneten Phase eine Divergenz der dielektrischen Konstante auf groBen 
Langenskalen auftritt. 

Fordert man im noch nicht renormierten Gitter die Giiltigkeit der ublichcn 
elektrodynamischen Beziehung 

V 4 x(r) = (3.1) 

so gilt offenbar £o = e(l) = in/ J. Mit Hilfe dieser Relation kann man eine 
renormierte Kopplungskonstante J(r) fiir beliebige Langenskalen definieren 

J(r) := iL. (3.2) 
e(r) 

Eine Multipolentwicklung || der Verzerrungsfunktion x(r) gibt AufschluB iiber 
die Struktur der dielektrischen Konstanten. 



3.1 Multipolentwicklung 

Mit der am Ende von Abschnitt 2.2 geschilderten Vorgehensweise kann man alle 
Versetzungen im System bei tiefen Temperaturen und kleiner Unordnungsstarke 
eindeutig zu Dipolen paaren. 

Fiir eine beliebige Lange £ konnen nun die Versetzungen in zwei Gruppen 
unterteilt werden. Die erste Gruppe besteht aus den Versetzungen, die zu Dipolen 
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mit einer Separation r < ( gebunden sind, die zweite aus denen, die grofiere 
Dipole bilden. Indiziert man die erste Gruppe mit T und bezeichnet die Versetzungen, 
aus denen der Dipol T besteht mit r^ und rr 2 , so liegt dessen Dipolzentrum 
bei 

r r := - (r Fl + r r J . (3.3) 

Die Quellenfunktion 77 (r) = eyVi Yl a b a ,j<>(r — r a ) unterteilt sich somit wie 
folgt: 



r?(r) =7 ?ex (r) + ^7 ?r (r) (3-4) 



wobei 



a 

lr(r) = e y Vj (b ri ,j 5(r - r r J + &r 2 ,j S(r - r r J) 
= eyVi br lt j (S(r - r r J - 5(r - r r J) . 



(3.5) 



Hierbei ist eine Summation iiber alle Versetzungen, die nicht zu Dipolcn 
mit einer Grofie r < C gebunden sind, also der zweiten Gruppe angehoren. Zu 
beachten ist, dafi die Anzahl der Versetzungen, die in rj ex {r) eingehen, von dem 
Wert von £ abhangt. Falls Q so grofi gewahlt ist, dafi alle Versetzungen Dipole mit 



r < C bilden, so ist rj ex (r) = 0. Nach (2.39) kann man die Verzerrungsfunktion 



nun ebenso wie die Quellenfunktion aufteilen x( r ) — Xex{r) + SrXr(r)- Der 
Beitrag des Dipols T zur Gesamtverzerrungsfunktion ist 

Xr(r) = J ( dVV(r')9(r - r'). (3.6) 



Im folgenden wird der Beitrag des Dipols Y zur Verzerrungsfunktion an einem 
gentigend weit entfernten Ort r mittels einer Multipolentwicklung abgeschatzt 
j9|. Diese wird fur jeden Dipol durchgefiihrt und die Beitrage der Entwicklungen 
schlicfilich summiert. Wie im Anhang A. 2 gezeigt wird entspricht dieses Vorgehen 
dem der Elektrodynamik. 

Zunachst ist es zweckmafiig, zu einer Darstellung mit inneren Koordinaten p' r ' 
des Dipols V iiberzugehen. 

Abbildung 3.1 verdeutlicht die vorliegende Situation. In inneren Koordinaten 
liegen die Versetzungen des Dipols T bei p\ = rr 1 — rr und p 2 = rr 2 — rr. 



In inneren Koordinaten wird aus (3.6) 

Xr(r) = J J dV r ^r(p (r) +r r ).g(r-r r -p (r) ). (3.8) 

Definiert man schlicfilich die Quellenfunktion mit ?7r(p < - r - ) ) := Vr(p^ +rr) um, 
so erhalt man 



(r) = J / d 2 p( r ),y r (p (r) )5(r-r r -p (r) ). (3.9) 
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4 

,<% 

i(r) 



i(r) 



bi 



Abbildung 3.1: Ein Dipol mit Zentrum bei it lieferi einen Beitrag zur 
Verzerrungsfunktion am Punkt r. Die innere Koordinate p^ r ' wird vom Zentrum aus 
gemessen. 



Da p( r ) eine interne Koordinate des Dipols ist, gilt in geniigend groBem Abstand 
vom Dipol |p( r )| <C |r — rp|. Somit kann obige Bczichung nach entwickelt 
werden. 



Xr(r)=J J dV r) ^r(p (r) )5(r-r r ) 
+ ij|dV r) ryr(p (r) )p (r) 4 P 



(r) d 2 g{r - r r ) 
3 dxndxrj 



(3.10) 



Hicrbci handelt es sich um eine ublichc Multipolentwicklung in zwei Dimensionen. 
Zu beachten ist, da8 man die Ablcitungen der Green- Funktion nach Komponcntcn 
von p( r ) durch Ableitungen nach rr ersetzen kann. Ein Term erster Ordnung 
der Form 



(3.11) 



tritt hier nicht auf. Einsetzen der Quellenfunktion fur einen Dipol zeigt sofort, 
daB er null ist. Wie bereits in Abschnitt 2.2.1 angemerkt, muB hier vcrglichen mit 
der Elektrodynamik doppelt so oft abgeleitet werden. Aus diesem Grund hat der 
Dipoltcrm obiger Multipolentwicklung die Form des klassischcn Quadrupolterms. 



Der Gesamtbeitrag aller Dipole mit r < ( zur Verzerrungsfunktion am Ort r ist 
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damit 

x ^(r)=^ Xr (r) 



j^(/dV r ), 7r( p(n )5(r _ rr) 

r •* 



(3.12) 



wobei 



^(r') =£ I cfp^fjr{p^)6(r'-rr), (3.13) 
r J 

C° lp (r') = £ / rf 2 P (r) ^r(p (r) ) P (T \ P (r) , 5(r' - r r ). (3.14) 
r ^ 

Mdchte man nun eine makroskopische Verzerrungsfunktion erhalten, so ist iiber 
Flachen AA 3> C 2 zu niitteln. Danach treten Dipole mit einer Separation r < C 
nur noch als polarisierbares Dielektrikum auf. Da eine Mittelung einer Grofie 
B(r) iiber eine Flache AA allgemein die Form (B(r)) = tq- f A 4 d 2 £ _B(r + £) 



bcsitzt, finden wir exemplarisch fur den ersten Term von (3.1 



= ^ J AA d 2 Z J dVy^(r"+05(r-r") 



(3.15) 



= J J rfV(^(r")) 9 (r-r"). 

Die hier auftretende Mittelung von ( |3.13| ) ist 

(V mp (r')) = 4rE / ^V r) ^r(p (r) ) / d 2 ^(r' + £-r r ). (3.16) 

Die (5-Funktion gibt fur einen Dipol T nur dann einen Beitrag, wenn dessen 
Zentrum rr innerhalb der Flache AA um r' liegt. Da AA ^> £ gilt, ist die 
Zahl dcr Dipole, die vollstandig in AA liegen wesentlich grofier als die Zahl der 
Dipole, die nur zur Halfte, das heifit mit nur einer Versetzung in dieser Flache 
liegen. 

Die bisher betrachtete Mittelung pro Flacheneinheit ist geeignet fur makroskopische 
Felder. Die Grofie rj Dlp (r') lafit sich aber auch pro Dipol mitteln. Ist AiV die 
mittlere Zahl von Dipolen mit r < £ in einer Flache der GroBe AA, so hat eine 
solche Mittelung hat die Form 



(B(r))AN = i/ A / 2fB(r+$) ' (3 ' 17) 
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Mit der makroskopischc Dipoldichte n£(r') bei r' erhalt man 
(r, Di P(r')) = (r, Di P(r')) AN n<(r'). 



(3.18) 



Schreibt man die makroskop. Verzerrungsfunktion x( r ) — Xex( T ) + [(x( r ))]o m 
der Form 



(r) w J / d 2 r'rf(v')g{v -r') + J dVcf,(r') 



dx'idx' 



(3.19) 



so ergibt sich fur die makroskopische Quellenfunktion nach einer zusatzlichen 
Mittelung iiber die Unordnung [...]d entsprechend (3.4), ( |3. 13 ) und ( [3. IS ) 



77 c (r) 



[ d 2 r'fjr{r')5(r-r r )\ n c (r) 
v r J I AN 



+ Vex(r). 



(3.20) 



Geht man fur den zweiten (3.12) analog zum ersten Term vor, erhalt man 
den Tensor CV, auf makroskopischem Niveau, der der elektrischen Polarisation 
entspricht 



4« = o 



J2 f cfr'firir'Wix'jSir-rr)) n c (r) 
\ r J I an 



(3.21) 



Partielle Integration und Ableiten liefern aus ( |3.19D 

V 4 X (r) = J V Hr) + J J d2r 'dJ^T C WW r r ') 



(3.22) 



J?7 C (r) + J 



dxidxj 13 



CUv). 



In der Elektrodynamik ist die Polarisation eine Funktion des elektrischen Felds 
und damit eine Funktion der Ableitung des elektrostatischen Potentials. Analog 
dazu ist der Tensor Cy- hier eine Funktion der zweifachen Ableitung der Verzerrungsfunktion 



Cj.(r) = £atfw(0 



;,/ 



d 2 x( r ) 
dx k dxi 



klrnn 



d 2 x(r) d 2 x(r) 

dx k dxi dx m dx n 



(3.23) 



Beriicksichtigt man nur den ersten Term und setzt ein isotropes Medium oder 
ein Dreiecksgitter an, so gilt in symmetrisierter Form 



d 2 x(r) 
d 2 x k 



(3.24) 



Diese spezielle Form von dijuiC) beruht auf der Tatsachc, dafi dieser Tensor 
- wie schon der Tensor Kijki aus ( |2.E| ) - den Gitterinvarianzen geniigen mufi. 
Eingesetzt mit (3.22) ergibt sich jetzt 



e(C)V 4 x(r)-47r^(r) 

47T 



= yV 4 x(r) - 47^(0 + C 2 (Q) V 4 X (r). 



(3.25) 
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Wir haben nun einen Ausdruck fur die dielektrische Konstante gefunden, namlich 
e(r)=e -^(C 1 (r) + C 2 {r)). (3.26) 

Ein Vergleich mit der iiblichen Formulicrung der Elektrodynamik zeigt, dafi 
— C\ — C 2 somit einer elektrischen Suszcptibilitat k cntspricht, die sich durch 
eine Dipolwahrscheinlichkcit P(r, 8) und eine Dipolpolarisierbarkeit a(r, 8) fiir 
Dipole mit Separation r und Winkel 8 wiefolgt ausdriicken lafit 

«(r) = -Ci - C 2 = / d 2 r'[P{r',8)a(r',8)] D . (3.27) 
Jo 

Ziel der folgenden Abschnitte ist es die relevantcn GroBcn zu bestimmen, so 
daB mittels eines expliziten Ausdruckes fiir die dielektrische Konstante e(r) der 
RcnormicrungsfluB der Kopplungskonstanten J(r) hergeleitet werden kann. 

3.2 Dipol-Wahrscheinlichkeit 

Im folgenden soli die Wahrschcinlichkeit P(r, 8) bestimmt werden, einen Versetzungsdipol 
mit Burgersvektoren b bei r a und -b bei zu finden. Der Winkel 6 ist wiederum 
der Winkel zwischen Abstandsvektor und Burgers- Vektoren, ferner ist r = 
\r a — 173 1 die DipolgroBe. 

Fiir das weitere Vorgehen beschranken wir uns auf die Untersuchung des Quadratgitters. 
Die meisten Argumentc konnen unverandert fiir das Dreiecksgitter ubernommen 
werden, notwendige Modifikationen werden in Kapitel 4 untersucht. Offcnbar 
cxisticren vier mogliche Orientierungcn fiir einen Dipol mit Burgers- Vektoren 
bei r a und r^, die in Abbildung 3.2 dargestellt sind. 




Abbildung 3.2: Schematische Darstellung der moglichen Orientierungen eines Dipols 
auf dem Quadratgitter 

Die mit 1 und 3 gekennzeichnetcn Einstellungcn haben beide die elastische 
Encrgie 

E = ^-lnr- -^-cos 2 8 + 2E c , (3.28) 

4-7T 47T 
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und eincn Bcitrag aus der Wechsclwirkung mit der Unordnung von V, beziehungsweise 
von —V. Fur die Einstellungen 2 und 4 ist 



E = E 



J , J ■ 2 

= - — In r — - — sm 

47r An 



2E C , 



(3.29) 



wahrcnd die Unordnungs- Wechsclwirkung hier mit V oder — V angesetzt wird. 
Bcriicksichtigt man die Moglichkeit, dafi kein Dipolpaar gebildet wird (Gewicht 
= 1, da keine Energie), so lautet die lokale Zustandssummme 



Z = 1 + w + + w_ + w + + w_ 
wobei die Gewichte durch 



w± 



-P(E±V) und 



W± 



-(3(E±V) 



(3.30) 



(3.31) 



gegeben sind. Somit ist die normierte Wahrscheinlichkeit ein orientiertes Paar 
mit dem Vektor 9 zwischen Burgers- Vektoren und Abstandsvektor zu linden 
(das sind die Einstellungen 1 und 3 aus Abbildung 3.2) 



P(r,l 



1 + W+ + W- + W+ + W- 



(3.32) 



D 



Da die Berechnung dieses Ausdrucks recht aufwendig ist, vergleichen wir inn mit 

der Wahrscheinlichkeit p(r, 8), die gegeben ware, gabe es nur die Orientierungsmoglichkeiten 

1 und 3. Dann ware die Wahrscheinlichkeit p 



U>4 



W- 



D 



Offcnbar ergibt sich direkt aus d3~32T ) und ( |3~35| ) 

P(r,6) <p(r,6). 



(3.33) 



(3.34) 



Da die Gesamtwahrscheinlichkcit Dipole mit Abstand r zu linden durch zusatzliche 
Einstellmoglichkeiten auf kcincn Fall absinken kann, mu8 auch die winkclgemittelte 
Dipolwahrscheinlichkeit P(r) mit vier erlaubten Burgers- Vektorrichtungen grofier 
sein, als die winkelgemittelte Wahrscheinlichkeit p(r) fiir nur zwei erlaubte 
Vektorrichtungen. Es ist 



271 



P(r) 
p{r) - 2 



d6 
2tt' 

2tt 



2tt 



d6 
2^ 



P{r,6) 



-[2P(r,9) + 2P(r,e + 7r/2)} =4 
dO , , 

Mit folgt hieraus die Ungleichung 

1 f 2n d9 . ns f 27r d9 nl m f 27r d9 . 
» / 5-p(r,fl)< / ^P(r,9)< / — p(r, 

2 J Q 2ir J 2tt J 2tt 

Da aber im folgenden Vorfaktoren nicht wesentlich sind (auch in || wurden nur 
die Hauptabhangigkeiten betrachtet), konnen wir somit als Naherung 



(3.35) 
(3.36) 

(3.37) 



P(r) (x p{r) 



dB 
2^ 



P(r,9) 



(3.38) 



ansetzen. 
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3.2.1 Wahrscheinlichkeit bei T = 



Setzt man die Gewichte ein, so liefert ( 3.38 ) 



l r°° dv _v e -e(E+v) + e -fj(E-v) 

P{r, 0) = -2 .__ e 2A 



(3.39) 

dV _ v* e -^E-v) 

: e 2A 2 



V / 27A 2 ~ " 1 + e-^f^-v) ' 

Die Summanden u> + konnen, da sie einen postiven Beitrag aus der Unordnungswechselwirkung 
haben, bei tiefen Tcmperaturcn allgcmcin vernachlassigt werden. Fur T = gilt 
insbesondere 



e -P(E-v) ( o falls E-V >0 



1 + e -^ E - v ) \ 1 falls E-V <0 
Somit hat man fur groBe r (lnr > <r und In 7' 3> Ec/J) 

pir, 0) = I - e 2^ 
Je V2^A* 



(3.40) 



(3.41) 



In r — cos 2 



wie in ( p. 73} ) aus Abschnitt 2.4 bereits hergeleitet worden ist. Da wir Flufiglcichungen 
aber nicht nur bei T = berechnen mochten, ist eine Untersuchung auch bei 
hoheren Temperaturen notwendig. 



3.2.2 Wahrscheinlichkeit bei T < T* 



Fiir T > gilt ( J3.40| ) nicht. Hier liefert auch der Fall E - V > einen Beitrag 
zur Integration 



rE dV -J^ -B(E-V) f°° dV — 



(3.42) 



= p(r, 0)+p T =o(r, 0), 
der zu berechnen ist. Mittels Substitution erhalt man fiir diesen Anteil 

p(r,0) = 4= e^ £+ ^ 2A2 [** e"' 2 *. (3.43) 

V 71 " I f3A 2 -E 

%/2A 

Allgemein erfiillt im Falle grofier Werte a; die Funktion $(x) := f^e^di 
die Relation 



i-*w-^s«^(i-s + -]. 
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die sich im Fall T < T* anwenden lafit, da (3 hinrcichcnd groB ist: 



P(r,9) 



(3A 
V2 



f3A 2 -E 
V2A 



e ±f3 2 A 2 -/3E 



2^ 



V2A 

[3A 2 -e' 



V2A _ 

Ja 2 



0E 



(3.45) 



E 



(3A 2 -E 



PT=o(r, 0), 



dabei wird der zweite Term der zweiten Zeile aufgrund der niedrigen Temperatur 
vernachlassigt. Mit ( 3.41| ) ergibt sich fiir 



f3A 2 

PM) = /3 A 2 _ E PT=o{r, 9) 



(3.46) 



Diese Formel ist offenbar eine Verallgemeinerung der Aussage fiir T — 0. Es 
zeigt sich, dafi im gesamten neu zu betrachtenden Bereich T < T* die fiir 
verschwindende Temperatur gefundenen Abhangigkeiten dominieren. Die Dipolwahrscheinlichkeit 
P(r) ist in dicscm Bereich somit 



, dO _ i _ 
P(r) tx 2 / — r 4 "«e* 
' 2tt 



■ cos" 



(3.47) 



denn es gilt e x cos = e a e a cos 29 , und die modifizierte Besselfunktion ist durch 
I (x) = £ e ^co s2 e denniert . 



3.2.3 Wahrscheinlichkeit bei T > T* 



Analog zum obigen Fall gilt auch im hohen Temperaturbereich (3.43). Hier 
strebt allerdings &t= — > 0, wahrend die untere Integrationsgrenze fur groBe 

Liingenskalen unbeschrankt ins Negative wachst. Wegen J_ e~ x dx = ^ gilt 
im diesem Fall 

p(r,6) oce^ A2 -^. (3.48) 

Da dieser Term den Beitrag von pr=o bei hohen Temperaturen iiberwiegt, kann 
letzterer vernachlassigt werden. Setzt man die bekannten Werte fiir E und A 
(aus (|2~70| )) ein, so resultiert abschliefiend fiir p(r, 0) 

-4h e T^ cos2f> fiir T<T* 

p(r,9) oc { (3.49) 
-44-( 1 -#) e 4ir( 1 -#) c ° s2e fiir T>T* 
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wobei die Bedingungen In r ^> a und In r ^ ^f- erfiillt sein miissen. 

Fur die winkelgemittelte Wahrscheinlichkeit P(r) einen Dipol mit Separation r 
zu linden, erhalt man daher 

f 2r-^e^/ (^-) fur T < T* 

P(r) oc I (3.50) 
[ 2 r -T^r( 1 -ir)e^r( 1 -irh (l - #)) fur T > T* 

Beide Proportionalitaten der Wahrscheinlichkeitsdichte liefern fur T = T* = 4p 
dasselbe Resultat, sie sind somit konsistent. Durch den Ansatz ( 3.32| ) fur eine 



normierte Wahrscheinlichkeit ergibt sich bei hohen Temperaturen zwar eine 
Boltzmann-verteilte Wahrscheinlichkeit, die || im gesamten Temperaturbereich 
angesetzt hat. Bei niedrigen Temperaturen T < T* , wo die Boltzmann-Verteilung 
zu Widerspruchen in Nelsons Arbeit fuhrt, liefert obiger Ansatz eine Fermi- 
verteilte Wahrscheinlichkeit. Aufgrund der Beriicksichtigung des fermionischen 
Charakters der Versetzungsdipole treten bei tiefen Temperaturen nun im Rahmen 
des Gultigkeitsbereichs obiger Formeln keine sehr grofien Dipolwahrscheinlichkeiten 
-P(r) 3> 1 mehr auf. Daher ist es gerechtfertig im folgenden dielektrischen 
Formalismus die Versetzungspaare als nicht-wechselwirkende Fermionen zu behandeln. 



3.3 Polarisierbarkeit 

Mit obiger Wahrscheinlichkeitsdichte ist es moglich, die elektrische Suszeptibilitat 
zu berechnen. Dazu wird der Dipoltensor Qu (r, 9) eines Dipols mit Abstand r 



und Winkel 8 benotigt, der nach (3.21) die Form 



Q ij (r,6) = ^ J d 2 r-q rfi {v)xiX 3 

(3.51) 

= -^£ik£jl{bk£lmXm + t>iek m x m ) 

annimmt. Hierbei ist rj r £ (r) die Qucllcnfunktion des Dipols. Die Wechselwirkung 
dieses Dipols mit einem von aufien angelegten homogenen Feld entspricht formal 
der Quadrupolwechselwirkung der zweidimensionalen Elektrodynamik, wobei 
aufgrund der Definition von Qij kein Vorfaktor auftritt. Diese Wechselwirkung 
ist 

W = ^Q ij E ij . (3.52) 



Der Tensor Eij symbolisiert die zweifache Ableitung der Verzerrungsfunktion 
Eij := ViVjX. Die Verzerrung, die \ darstellt, wird durch eine von aufien an 
das System angelegte Konfiguration hervorgerufen. Analog zu Elektrodynamik, 
wo man Konfigurationen betrachtet, die zu einem homogenen, konstanten E- 
Feld fuhren, sollen die Konfigurationen hier so gewahlt sein, dafi der Tensor 
konstant ist. 
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Auch hier geniigt es, sich bei der Untersuchung auf den Fall mit zwei erlaubten 
Burgers- Vcktorrichtungen zu beschranken, da die hier untcrsuchtc Polarisicrbarkeitsdichtc 
und damit die Suszeptibilitat durch die Einfiihrung zusatzlicher Ausrichtungen 



nicht abnehmen darf. Es gilt somit eine zu (3.37) analoge Relation fur diese 
Grofien. 

Die Polarisationsdichte qij{r,9) fur einen Versetzungsdipol mit Burgers- Vektor 
b bei r Q und mit -b bei rg ist somit 



Qij(r, Q) 



w+ — W- 
+ w + + w 



-Qij(r,0) 



D 



mit 



y± = exp [-/3(E ±V±W)] 



(3.53) 



(3.54) 



Mittels dieser Polarisationsdichte lafit sich formal eine Grofie Ayjy definieren 



aus deren Komponenten dann mit (3.24) die Polarisierbarkeitsdichte der Versetzungspaare 
bei Abstand r ermittelt werden kann. 



Aijki (r) 



1 



2tt 



d9 


~dqij(r, 6)' 




2^ 


[ dE M 


D 



(3.55) 



Diese Definition entspricht dem Vorgehen in der Elektrodynamik, wo ein ebenso 
angesetzter Tensor die Bercchnung der Polarisierbarkeit erlaubt. Diese Analogic 
wird ausfuhrlich in Anhang C.l gezeigt. 

Der Vorfaktor | ist notwendig aufgrund der Tatsache, dafi in ( 3.53| ) eine Orientierung 
nur modulo 180 Grad ausgezeichnet ist, der Faktor 2, da es auf dem Quadratgitter 
neben den beiden in ( 3.53} ) beriicksichtigten Burgers- Vektor - Orientierungen mit 
Winkel 8 noch zwei weitere mit Winkel 9 + ^ gibt, die denselben Beitrag zur 
Polarisationsdichte liefern. 



3.3.1 Polarisierbarkeit ohne Unordnung 



Zunachst werten wir (3.55) fur den Fall ohne Unordnung aus. Hier erhalt man 
nach der Ableitung 



Aijki (r) 



d9 
2^ 



PQijQkl[(w+ + + w + + w-) - (w+ - w^) 2 } 



(1 + W+ + W-) 2 



2-7T 

— (3QijQkl 



w + + w_ 
1 + w + + u;_ 



_E=0 



J ^PQvQ ki2 p( r >< 



(3.56) 



dabei wurde benutzt, dafi ohne Unordnung und ohne aufieres Feld w + — w- = 
gilt. Wie in Anhang C.2 hergeleitet wird, gilt folgende niitzliche Relation 
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^ de 



QiiQki e x cos 



/o 2- lJ 



1 X 



(3.57) 



wobei 



<3,M = ^/„(f) 

<J»(r>~±^.(§) 



(3.58) 



Jo und I\ sind modifizierte Besselfunktionen, deren Definition durch 

f 2-tt fja 



M*) := / 



£ cos 28 



(3.59) 



M^) : = / ^-cos20e a 
Jo 27r 

gegeben ist. Definiert man nun die GroBen Ai(r) und ^(r) wie folgt 

d{r)= [ d 2 r'Ai(r') miti = (l,2), (3.60) 



o 



so resultiert aufgrund von (3.24) und der speziellen Form der Dipol-Wahrschein- 
l ichkc it ohne Unordnung (p(r,0) — r^Th' e^r cos s ), die die Benutzung von 
( |3.57i ) ermoglicht: 

— ^4^(r) A X y X y{v^) A X yy X (T^ Ay XX y{v^j Ay X y X \T) 



2 1 j_ _j_ r 2 (J 

■ r 4 " T e 8 " T — 7 



T 4 2 u \8ttT 



2A 2 (r) = A xxxx {r) + A xxyy (r) - A\(r) 



t 4xT ^4, i ( t ( J \ t ( J \ T f J 



(3.61) 



e^r 2 I —)-h[ — )- Io 



2T \ \8ttT J \8ttT J \8ttT 



1 ■]_ _J_ o T I J 

= r 4,rT e S7vT r z / x 



2T V 8ttT 



40 



3.3. POLARISIERBARKEIT 



Mit Gleichungen ( gjg ) und fiTTft zcigt sich nun, da8 die Polarisierbarkeitsdichtc 
A(r) folgenden Wert annimmt 



Mr) 



(A 1 (r)+A 2 (r)) 



T" ■]_ _J_ . 



8ttT 



V— 

2 V8ttT 



(3.62) 



Diese Polarsierbarkeitsdichte stimmt mit den von Kosterlitz und Thouless Q 
angegebenen Werten fur C\ und C% iiberein. 

Fur die Polarisierbarkeit a(r) eines Versetzungspaares folgt unter Beachtung 
der Relation Q3.27D 



a(r) 



A(r) 



Jd9^ (1- ± cos 26 



f 4ttT e 4ttT 



7" 4ttT g 4-rrT 



S 2 



(3.63) 



2T 



< cos 20 > c 



Ohne Unordnung ist dieser Wert fur die Polarisierbarkeit bereits von Young |l3| 
auf eine etwas andere Weise hergeleitet worden. Das Symbol < ... > ang bedeutet 
eine Mittelung iiber die Winkcl. Der Wichtungsfaktor ist die Dipolwahrscheinlichkeit 
und damit vom Temperaturbereich abhangig. 



3.3.2 Polarisierbarkeit mit Unordnung 



Dieses bis hierher entwickelte Prinzip zur Behandlung der Polarisierbarkeit 
ohne Unordnung kann durch eine Erweiterung auch auf Falle mit Unordnung 
iibertragen werden. Hierzu ist eine Unordnungsmittelung in Formel ( |3.56} ) notwendig 
Der Tensor A^i (r) wird modifiziert zu 



A ijk i(r) 



2w dO 

7T~ PQijQkl 

Zn 



w + + w_ + Aw + W- 
(1 + w+ + W-) 2 



E=0. 



(3.64) 



:=7r(r,0) 



Es ist also gegeniiber dem bisher Berechneten zusatzlich der Ausdruck n(r, 9) 
auszuwerten. Beschrankt man sich dabei in Zahler und Nenner auf die fuhrenden 
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Terme, so folgt 
7r(r, 0) w 2 



-0(E-V) 



o \/2ttA 2 



1 + e -2^(B-V) 



-/3(B-^) 



o V2ttA 2 



V27A 2 



(3.65) 



Unter Benutzung der Methoden und Schrcibwcisc aus Abschnitt 3.2.2 laBt sich 
dies vereinfachen zu 

7r(r,fl) =2p(r,0) + ^e^ + ^ 2A2 / e^dt 

i/2A 



2p(r, 9) + e 



/3_B+i/3 2 A 2 



(3A 2 + E 
V2A 



$ (oo) 



(3.66) 



= 2p(r, 6») + 2 



E + f3A 2 



PT=o(r, 0). 



Im hohen Temp eratu rbereich j 1 > T* gilt p 3> Pt=o, daher kann man den 
zweiten Term in ( |3.66 ) vernachlassigen. Ein identisches Vorgehen wie im vo rherig en 
Abschnitt ist moglich, da auch hier die Form von p die Anwendung von ( 3.57 ) 
erlaubt. Wir erhalten fur die Polarisierbarkeitsdichte und die Polarisierbarkeit 



4TJ-J5 AT J — .15 
T 



Air) = T WttT' z e 32ttT 



UTJ-Ja \ 1 (\TJ-Jo 
I 32ttT 2 ) ~ 2 1 V 32ttT 2 



a(r) 



2T 



1 - - < cos 29 > 



ang> 



(3.67) 
(3.68) 



Fiir Temperaturen unterhalb der Grenze T* gibt Relation ( 3.46| ) 

E E 



(3 A 2 - E + E + [3 A 2 
T 

■PT=o(r, 0). 



PT=o{r,< 



T + T 



(3.69) 



Nach Anwendung des bekannten Verfahrens findet man auch hier die gesuchten 
GroBen. 



Insgesamt ergibt sich so fiir die Polarisierbarkeitsdichte A(r): 



T+T- 



L l 

■ f 47TCT g StTO 



fiir T <T* 



A(r) = 



r „— 
2T 



ATJ-J5 ATJ-Ja 



r- We W [J (^) - ij, (^)] fiir T>T* 

(3.70) 



42 



3.3. POLARISIERBARKEIT 



Analog ist die Polarisierbarkeit 



fl-i <cos20> fflB „^) fur T<T* 



T+T* V 2 yAJo aw ^ang<) 
n(r) = { ' (3.71) 

^(l-± <cos20> a „ s> ) fur T>T* 

Die Notationen der Winkclmittelung sollen den Unterschied der zu benutzenden 
Gewichte verdeutlichen. Es sind 

"■ an9< 3 ^ , (3.72) 

^ ■■■ ^ang> ■ c 

Diese Resultate sind konsistent bei T = T* . Die angestrebten Flufigleichungen 
bei Renormierung konnen jetzt hergeleitet werden, indem man die gewonnen 
Ergebnisse in (3.26) einsetzt und ( |3.2D beachtet. Das nachste Kapitel wird diese 
Flufigleichungen und ihre physikalischen Konsequenzen behandeln. 
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Kapitel 4 



Flufigleichungen und 
Phasendiagramme 

4.1 Ergebnisse fur das Quadratgitter 

In Kapitel 3 sind alle notwendigen Relationen hergeleitet worden, um das Schmelzen 
kristallincr Filmc mit quadratischem Gittcr zu analysicrcn. Das Skalenverhaltcn 
der dielektrischen Konstante und damit dcr Kopplungskonstante bei Renormierung 
kann aus den Ergebnissen gewonnen werden. Dieses Skalenverhalten la8t direkte 
Schliissc auf das gcsuchtc Phasendiagramm zu, welches die feste Phase mit quasi- 
langreichweitiger Ordnung beziiglich Translation von einer nicht-kristallinen 
Phase trennt, die diese Ordnung nicht zeigt. 



4.1.1 Herleitung der Flufigleichungen 



Zunachst ist durch ( 3.26 ) der Wert der renormierten dielektrischen Konstanten 

in Abhangigkeit der Langenskala gegeb en. Mit der Beziehung zwischen Suszeptibilitat 

zur Polarisierbarkeitsdichte aus ( 3.60| ) gilt daher 



e(r + dr) = e(r) + A-kA{t) 2irr dr. (4.1) 

Nun ist es moglich die Anderung der Kopplungskonstanten J(r) zu ermitteln, 
die resultiert, wenn man bei der Renormierung von einer Langenskala r zu einer 
Skala r+dr iibergeht, das heifit den Abschirmeffekt derjenigen Verse tzu ngspaaren 
zu beriicksichtigen, deren Separation zwischen r und r + dr liegt. ( |3.2| ) gibt mit 
Hilfe einer Entwicklung 

J(r + dr) 



e(r + dr) s(r) + 8it 2 A(r) r dr 

(4.2) 

J(r) - J(r) 2 A(r) 2nr dr + 0(A 2 ). 
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Hieraus folgt sofort die FluBglcichung der Kopplungskonstanten. Setzt man I := 
lnr, so ist 

^-=r^- = -J(r) 2 A(r)2nr 2 +0(A 2 ). (4.3) 
dl dr 

Wir fiihren nun den bereits in Kapitel 2 behandelten Begriff der Fugazitat y 
ein. Wegen ( 2.7l| ) ist der Ausgangswert dieser Fugazitat im nicht-renormierten 



System mit r = 1 gegeben durch 

E c (C + 1)J 

y = e ^ = e s.r ( (4.4) 

die Konstante C betragt nach Anhang A etwa 5. In renormierten Systemen 
wird die Fugazitat als y — e~ Fc l T aufgefafit, wobei Fc die freie Energie eines 
Versetzungskerns auf der entsprechenden Skala r ist. Da die Zahl der durch 
Renormierung verlorengegangenen moglichen Dipolplatze pro Platzpaar oc r 4 
skaliert und die Wahrscheinlichkeit einen Versetzungsdipol zu finden p(r, 9) 
betragt, ergibt sich fiir die renormierte Fugazitat 

y(r) 2 = < r 4 p(r,e) > . (4.5) 

Hierbei bezeichnet < ... > eine naher zu bestimmende Art der Winkelmittelung. 
Beriicksichtigt man fiir hohe Temperaturen T > T* , den bisher vernachlassigten 
Beitrag der Energie der Versetzungskerne Ec zur Dipolwahrscheinlichkeit, so 
modifiziert sich ( 3. 49] ) zu 



p(r,0)cxr-^( 1 -^)e^( 1 -^) cos29 e -^. (4.6) 



Dies zeigt, dafi (|f.5|) fiir kleine r der Forderung yo — e Ec ' T nur dann geniigcn 



kann, wenn die Winkelmittelung die bereits fiir die Polarisierbarkeit(vgl. (3.72)) 
benutzte Form hat: 

y(r) 2 =<r 4 p(r,d) > ang . (4.7) 

Dicsc Definition der renormierten Fugazitat stimmt mit den von Nelson, Halperin 
und Young genutzten Ansatzen O, |l^| iiberein. Fiir tiefe Temperaturen T < 
T* findet man yo = e~ Ec l T allerdings nicht, wenn man die Dipolwahrscheinlichkeit 
entsprechend ( 3.4£| ) ansetzt, welche aber auch nur fiir grofie Werte von r korrekt 



ist. Eine genauere Untersuchung fur kleine Langenskalen gibt auch hier das 
richtige Ergebnis. 

Schliefilich bleibt noch zu bemerken, dafi mittels dieser Definition die Fugazitat 
im relevanten Bereich a < a und T < T* klein bleibt. Hier gilt auf grofien 
Skalen \y\ <C 1. Wir fiihren jetzt noch eine niitzliche Variable r* ein, die im 
folgenden die expliziten Fallunterscheidungen zwischen hohem und niedrigem 
Temperaturbereich uberflussig macht 

2T 

t* = min(l, — ) (4.8) 
Ja 

Somit nimmt r* im Bereich T > T* den Wert 1 an, im Bereich T < T* den 
Wert t* = Insbesondere ist diese Definition konsistent bei T = T* . Fiir die 
renormierte Fugazitat erhalt man hicrmit auf grofien Skalen 

y(/) 2 =r 4 r-^( 1 -^ 1 ) . (4.9) 
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Die hier sowie im folgenden auftretenden Grofien sind immer renormiert. Mit 
dieser Fugazitat und ( |4.3} ) ergeben sich die gesuchten Renormierungs-Flufigleich- 
ungen fiir das quadratische zweidimensionale Gitter: 

2nd 2 



dJ 
~dl 



dy 
dl 



do 



T 



To 



-T(t*) 

(— 

\8ttT 



V exp 



/ Jt* 
\8^f 



Jar" 
AT 



2- 



Jt* 
8^T 



1 - 



Jar* 
~4T~ 



=0 



1 



1 

-2 h 



Jar* 
4T 
f Jt* 



\8ttT 



1 - 



Jar" 

It" 



(4.10) 

(4.11) 
(4.12) 



Bei diesen FluBgleichungen wurden lediglich die ersten beiden Ordnungen in 

y beriicksichtigt, was aufgrund obiger Ausfiihrungen gerechtfertigt ist. Eine 

Renormierung der Unordnungsstarke kann mit dem angewandten Formalismus 

nicht gefunden werden. Die von Tang || und Scheidl J(| gefundene Unordnungsrenormierung 

am XY-Modell miisste bei Verwendung der von ihnen genutzten Formalismcn 

auch bei diesem Problem auftretcn. 

Die FluBgleichungen stimmen fur Systeme ohnc Unordnung mit den Ergcbnisscn 
von Nelson, Halperin und Young iiberem Jl^, . Den Fall mit Unordnung hat 
Nelson Q lediglich am Dreiecksgitter untcrsucht, aber eine Anwendung der 
dort benutzten Methodik auf das Quadratgitter gabe eine Ubereinstimmung im 
gesamten Bereich T > T* . Bei tiefen Temperaturen wiirde sich mit Nelsons 
Methodik auch fiir das Quadratgitter ein Wiedereintritt in die nicht-kristalline 
Phase ergeben, den unsere Gleichungen nicht aufwciscn. 



4.1.2 Flufi und Phasendiagramm 



Wertet man die Flufigleichung der Fugazitat, so ist folgendes Verhalten offensichtlich 

(4.13) 



< falls T <T* und a < r^- 
dy j oderT>T*undCT< - ^ L 



> sonst 

Wie bereits in Abschnitt 2.3.1 erlautert, kann in einem Bereich, in dem y auf 
grofien Skalen immer weiter ansteigt, keine geordnete Phase vorliegen. Aufgrund 
der Proportionalitat zwischen y 2 und der Dipolwahrschcinlichkeit, wird diese 
dann fiir r — > oo sehr grofi, was gleichbedeutend mit der Existenz freier Versetzungen 
ist. Eine Bedingung fiir die feste Phase ist also ^ < auf grofien Langenskalen. 

Mochte man den durch ( gig ) bis ( p^ ) erzeugten Flufi analog zu Abbildung 



2.3 veranschaulichen, so ist es zunachst notwendig 

dy dy dl dJ f T ^ 2 



dK- 1 dldJdK 
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mit K := ^ zu berechnen. Im Bereich 1 > § ist keine Veranderung des von 
Nelson gewonnen Flufibildcs (vgl. Abbildung 2.3) zu erwarten, da die Flufigleichungen 
im Bereich hoher Temperatur unverandert sind. Intcressant hingegen ist das 
Verhalten im tiefen Temperaturbereich. 

Numerische Integration von ( 4.14| ) liefert fiir einen festen Wert von a den in 
Abbildung 4.1 dargestellten Renormierungsflufi. 





Abbildung 4.1: Renormierungsflufi des Quadratgitters bet a — 



Der Startpunkt des Flusses befindet sich immer auf der gestrichelten Linie yo- 
Die dick eingezeichnete Linie ist die Separatrix, die einen schraffierten Bereich 
abtrennt. Sie miindet bei KT in die Linie y = 0. Der Schnittpunkt der Separatrix 
mit der Linie yo liefert einen Wert K~ x . Fiir Anfangswerte mit K" 1 < K~ l , also 
bei geniigend tiefer Temperatur, liegt der Startpunkt im schraffierten Bereich, 
wo der Flufi immer bei y = und einem endlichen Wert von K endet. Hier 
befindet man sich somit in der festen Phase. Liegt der Startpunkt aufierhalb 
des schraffierten Bereichs, wachst die Fugazitat immer starker an, man befindet 
sich in der nicht-kristallincn Phase. 

Somit ist T c (a) = K~ 1 J, die reale Phaseniibergangs- Temperatur, sie fallt fiir 

Ec — * oo mit T + zusammen. Analysiert man fiir alle Werte von a den Renormierungsflufi, 

so ergibt sich das Phasendiagramm fiir das System, dargestellt in Abbildung 4.2. 

Die durchgezogene Linie stellt die renormierte Phasengrenze T + (a) dar, wah- 
rend die reale Phasengrenze T c (d) durch die gestrichelte Linie wiedergegeben ist. 
Im Gegensatz zum Phasendiagramm von Nelson (dim gestrichelte Linie) findct 
man hier keinen Wiedereintritt in die ungeordnete Phase. Daher steht dieses 
Ergebnis auch nicht im Widerspruch zu der Arbeit von Ozeki und Nishimori 

HI- 
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-J- T 

16tt 



Abbildung 4.2: Phasendiagramm fiir das Quadratgitter. Es zeigt die Abhangigkeit 
zwischen kritischer Temperatur und Unordnungsstarke 



4.2 Ergebnisse fur das Dreiecksgitter 

Der wesentliche Unterschied zwischen Quadrat- und Dreiecksgitter, der dazu 
fiihrt, daB die Ergebnisse des Quadratgitters nicht einfach auf das Dreiecksgitter 
iibertragbar sind, ist die Tatsache, daB auf dem Dreiecksgitter auch drei Burgers- 
Vektoren so kombiniert werden konnen, daB sie in der Summe ladungsneutral 
sind. Wir werden diesen Unterschied mit Hilfe des von Young Jl3| benutzten 
Formalismus behandeln. 



4.2.1 Modifikationen fiir das Dreiecksgitter 

Eine solche ladungsneutrale Kombination von drei Burgers- Vektoren, im folgendcn 
auch Tripol genannt, zeigt Abbildung 4.3. Der Abstand r' sei der kleinste 
Abstand zweier Burgers- Vektoren. Uberschreitet die Renormierung nun die Lange 
r' , so werden die Vektoren 1 und 2 zu einem Burgers- Vektor zusammengefafit. 
Auf groBeren Langenskalen liegt somit effektiv ein Dipol vor. Tripole konnen 
daher in die Berechnung einbezogen werden, in dem man erlaubt, daB einer 
der beiden Burgers- Vektoren jedes Versetzungspaares mit Separation r aus zwei 
Burgers- Vektoren mit einem Abstand r' < r zusammengesezt ist. 

Offcnbar andert sich hiermit die Dipolwahrscheinlichkeit, was Neuuntersuchungen 
crforderlich macht. 

Zunachst ist es notwendig die elastische Energie ohne Unordnung dicser Konfiguration 
zu ermitteln. Wie auch fiir Dipolc auf dem Dreiecksgitter, hat ein Tripol bei 
festen Postionen sechs Orientierungsmoglichkeiten fiir die Burgers- Vektoren. 
Abbildung 4.3 zeigt eine davon. Fiir den Fall r 3> r' findet man aufgrund der 
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Abbildung 4.3: Darstellung eines Tripols. Fafit man die Vektoren 1 und 2 zusammen, 
so erhalt man effektiv einen Dipol mit Separation r und Wmkel 9 



Relation ( |2.45| ) 



E Tr t = 2 ^- lnr + — lnr' + 3E C + ^-cos0 cos(6> + ^) 
8n 8n 47r 3 

J 2-7T J , 2ti\ 

H cost* cos(# ) H cos0 cos(0 H ). 

An 6 in 3 

Einfachc Umformungen mit 7 := ip — 6 ergeben 

E Tri = — In r + — In r' + 3E C - — cos 2 6 
An 8n An 



J 



J 



— COS- 7 + T7j— 

47r Ion 



(4.15) 



(4.16) 



— In r' cos 2 7 + Ec , 

8n in 



wobei E Dip die bereits in (|2.7l|) angegebene elastische Energie eines Versetzungsdipols 
mit Abstand r ist. 

AuBer der elastischen Energie mu8 auch die Varianz der Wechselwirkung des 
Tripols mit der Unordnung bereclmet werden. Diese Varianz ist 



tf Tri =[{V{vi) + V{r 2 )+V{T3)?]^ 



(4.17) 



Bei Berechnung dieses Ausdrucks ist zu beriicksichtigen, dafi die Burgers- Vek- 
toren hier nicht, wie bei Dipolen antiparallel zueinanderstehen, sondern der 
Winkel zwischen ihnen jeweils ^ betragt. Ist A 2 Dip = A 2 (r,8) die Varianz der 
Unordnungswechselwirkung eines Dipols entsprechend ( 2.70D , so gilt wegen der 
Burgers- Vektor-Stellung 



A 2 



A 2 Dlp + [V{v 1 ) 2 ] D + 2[V{v 2 )V{v 1 )] D . 



(4.18) 
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Aufgrund der Berechnungen aus Anhang B erhalt man 

At„ = &Dip + ^ lnr ' + ^ C0S <t> COS(0 + y ). (4.19) 



Formt man dies analog zu ( |4.16| ) um, so resultiert abschlicBcnd 

— lnr' - — cos 2 7 + — . (4.20) 



J 2 a , . J 2 a 2 J 2<J 
i '--^ + T6^ lnr -^T C0S 7+ 32^' 



Die modifizierte Dipolwahrscheinlichkeit setzt sich nun aus dem reinen Dipolanteil 



und dem Tripolanteil zusammen. Entsprechend (3.49) gilt 
p(r, 9) oc 



ex P (-iSfc) + ex P ("2§^) ffir T <T* 



exp (i^A^ p - + exp (|/3A|, ri - /3£ TH ) fiir T > T* 

(4.21) 

Setzt man alle bekannten GroBen ein und benutzt den Parameter t* , so ist 
p(r,6) = r~^( 1 -^)e^( 1 - d i^) c ° s2e (l + e^C 1 "^ 1 ) q(r)) . (4.22) 

Hierbei gibt g(r) die Zahl der Moglichkeiten an, dafi einer der beiden Burgers- 
Vektoren (daher Faktor 2) des Dipols aus zwei Vektoren kleineren Abstands 
kombinicrt ist: 

q(r) = 2 JJ J^d^dr'r'r''^ 1 '^^ e^ 1 '^-)^^ (4.23) 

Indem man erlaubt, daB einer der Vektoren, die den schon zusammengesetzten 
Dipolvektor bilden wiedcrum aus zwei Vektoren noch geringerer Separation 
gcbildet ist, kann dieser Betrachtung die von Young geforderte Selbstkonsistenz 
gegeben werden. Diese Art der Einfuhrung von Quadrupolen etc. fiihrt allcrdings 
lcdiglich zu Korrekturen hoherer Ordnung, die im folgenden nicht von Interesse 
sind und auf die daher verzichtet werden kann. 

Soil wie fiir das Quadratgitter die Relation 

y(r) 2 =<r 4 p(r,9) > ang (4.24) 
gelten, so ist die Fugazitat y fiir das Dreiecksgitter wie folgt zu definieren 

y(l) 2 ^rV^f 1 ^) (l + e-T^ri 1 -^) q(r)) . (4.25) 

Bei Untersuchung kleiner Langenskalen gilt auch bei dieser Definition der geforderte 
Anfangswert 

ya = e-^. (4.26) 

Somit sind effektive Dipolwahrscheinlichkeit und Fugazitat dem Dreiecksgitter 
angepafit worden. Die modifizierten Flufigleichungen folgen hieraus sofort. 
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4.2.2 Flufigleichungen fur das Dreiecksgitter 



Wie schon im Falle des Quadratgitters entwickeln wir auch hicr die Flufigleichungen 
bis zur zweiten Ordnung in y. Mit den Ableitungsregeln ergibt sich fur den 
Renormierungsflufi der Fugazitat 



<hi_ 
dl 



r dy 2 
2y dr 



(4.27) 



Arbeitet man mit der modifizierten Fugazitat, bleibt der Flufi der Kopplungskonstante 
in den betrachtctcn Ordnungen formal unvcrandert. Es erfolgt lediglich die 
Anderung des Vorfaktors, da in ( |3.55| ) der erste Faktor eine 3 statt einer 2 
sein mufi, befmdet man sich doch jetzt auf dem Dreiecksgitter, wo es insgesamt 
sechs Orientierungsmoglichkeiten fur einen Dipol mit festen Positionen gibt. Die 
Renormierungs-Flufigleichungen fur das Dreiecksgitter nehmen somit folgende 
Form an: 



dJ 



3nJ 2 



T + T{t* 
( Jt* 



2 ( Jt* 

— y ex P 



\8nT 



\8ttT 



1 - 



Jar" 



Jot* 
AT 
Jt* 
8^f 



Jar* 
AT 



(4.28) 



dy 
dl 



^= 2 



Jt* 
~8^T 



2tt exp 



Jar* 



— 

\16ttT 



AT 
1 - 



y 

Jar* 
AT 



( Jt* 



\8irT 



Jar* 
AT 



(4.29) 



d i=o 



(4.30) 



Wie bereits erwahnt ist mit dem verwendeten Formalismus keine Renormierung 
der Unordnungsstarke festzustellen, hierzu gelten die schon fur das Quadratgitter 
gemachten Bemerkungen. Wie schon die Flufigleichungen des Quadratgitters, 
stimmen die des Dreiecksgitters mit den von Nelson, Halperin und Young gefundenen 
Gleichungen fur Systeme ohne Unordnung iiberein. Zusatzlich ist die Ubereinstimmung 
mit Nelsons Resultaten fur ungeordnete Systeme Q im Bereich T > T* zu 
erwahnen. Ebenso wie im Fall des Quadratgitters aber zeigen obige Flufigleichungen 
auch fur das Dreiecksgitter nicht den von Nelson vorhergesagten Wiedereintritt 
in die ungeordnete Phase. 



4.2.3 Phasendiagramm fur das Dreiecksgitter 



Aufgrund des zweiten Summanden in der Flufigleichung der Fugazitat kann ^j- 
hier nicht so einfach abgeschatzt werden wie im Fall des Quadratgitters. Offcnbar 
bleibt (4.13) aber giiltig fur y — > 0, so dafi K+ 1 (K := ^) nicht verandert wird 
und keine Anderung der renormierten Phasengrenze T + (a) eintritt. 
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Dennoch ist es notwendig den Flufi gcnaucr zu untersuchen, so daB cr iiber 

numerische Intergration d ^ti veranschaulicht werden mu8. Geht man analog 

zum Quadratgitter vor, so erhalt man den in Abbildung 4.4 dargestellten Renormierungsflufi. 

Auch hier ergibt sich im Bereich K~ x > § keine Anderung gegenuber Nelsons 

Ergebnisscn. 




Abbildung 4.4: Renormierungsflufi des Dreiecksgitters bei a — 

Auch hier befindct man sich nur fur Startpunktc im schrafficrtcn, durch die 

dicker cingczeichnete Separatrix abgetrenntcn Bereich, in dcr kristallincn fcstcn 

Phase. Fur Startpunkte auBerhalb dieses Bereichs endet der FluB bei beliebig 

grofien Fugazitaten, also in der nicht-kristallincn Phase. Da allc rcalen Startpunkte 

auf der gestrichelten Linie y liegen, gibt der Schnittpunkt K~ x die reale Phaseniibergangs- 

Temperatur. 

Ein auffallcndcr Untcrschicd zum Renormierungsflufi des Quadratgitters ist, daB 
der FluB fur K^ 1 — > nicht im Unendlichen, sondern bei einem festen, nur von 
a abhangigen Wert, namlich bei 

^'-'-(sf-O^SibW (4 - 31) 

beginnt. Dies hat keine direkte physikalische Konsequenz, da der Wert fur a < 
— immer grofier als null ist und daher kein Wiedereintritt in die ungeordnete 
Phase moglich ist. Dennoch stellt sich die Frage, ob dieses Verhalten physikalisch 
ist, oder nur durch zu grobe Abschatzungcn in Abschnitt 4.2.1 hervorgerufen 
wird. 

Analysiert man hier den den RenormierungsfluB fur alle Werte von a, so erhalt 
man das in Abbildung 4.5 wiedergegebene Phasendiagramm des Dreiecksgitters. 

Wicderum stellt die durchgezogene Linie die renormierte Phasengrenze T+(a) 
dar, die gegenuber dem Quadratgitter unverandert bleibt. Lediglich die reale 
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Flussigkristall 




Abbildung 4.5: Phasendiagramm fur das Dreiecksgitter. Es zeigt die Abhdngigkeit 
zwischen kritischer Temperatur und Unordnungsstdrke 



Phasengrenze T C A (a) liegt bei einer tieferen Temperatur als die beim Quadratgitter 
T ( P(ct). Ein wesentlicher Unterschied zwischen den Phasendiagrammen des Quadrat- 
und des Dreiecksgitters kann also nicht festgestellt werden. In beiden Fallen 
erfolgt kein Wiedereintritt in die ungeordnete Phase. 



4.3 Entropieflufi 

Abschliefiend stellt sich die Frage, ob die gewonnen Ergebnisse zu einem nicht 
negativen und damit physikalischen Entropieflufi fuhren. Wie bereits in Abschnitt 
2.3.2 angemerkt, kann der Beitrag zur freien Energiedichte dFv , den Versetzungsdipolc 
mit einer Grofie zwischen r und r + dr leisten, durch 

»2tt 



1 f* 

dF v « -Trdr- J dd[\nZ] D (4.32) 



angenahert werden. Hierbei ist Z die lokale Zustandssumme, die fur das Quadratgitter 
durch ( fc.3C| ) gegeben ist. In diesem Fall ist somit 

1 f 27r 

dF v --Trdr- dd [ln(l + w+ + W- + w + + ui-)] D , (4.33) 
4 Jo 

da vier mogliche Dipo lcinst ellung en ex istieren. Deren elastische Energien E 



und E (gegeben durch ( 3.2§|) und ( 3.29 )) lassen sich nun in abstandsabhangige 



Anteile E r und E r , sowie in winkelabhangige Anteile E w und E w aufteilen, so 
daB E = E r + E w und E = E r + E w gilt. Auf geniigend grofien Langenskalen 
In r 3> ^f- erhalt man so 

E r = E r = — lnr (4.34) 

47T 
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und 



E w = -—cos 2 
4ir 



E w = -— sin 

47T 



(4.35) 



Ablciten von nach e~ pBr liefert 



, Er d{dF v ) 1 
S(e-^) Jmr 4 7(1 



■ w 4 



1 + w+ + u>_ + w+ 
Mit den Berechnungcn aus Abschnitt 3.2 erhalt man sofort 



D 



-fJE r 9{dFy) _ _7T 



—Trdr P(r), 
2 v h 



(4.36) 



(4.37) 



wobei P(r) die durch ( 3.5C| ) gegebene winkelgemittelte Wahrscheinlichkeit ist, 

einen Dipol mit Separation r bei vorgegebenen Versetzungspositionen zu findcn. 

Fiir das Dreiecksgitter gilt diese Relation auch, allerdings ist hier fur die Wahrscheinlichkeit 

P(r) der sich aus (4.22) ergebende Ausdruck einzusetzen. 



Im Fall T > T* = If folgt daraus 



dF v = ~-TrdrP(r), 



(4.38) 



da die P(r) hier linear in e l3Er sind. Fiir den FluB der freien Energiedichte 
ergibt sich 



dF v 
dl 



-Tr 2 P{r) 



(4.39) 



r 2 \8ttT \ AT 



Dieses Resultat entspricht, abgesehen von Vorfaktoren, den von Tang || und 
Scheidl Q gefundenen Ergebnissen fiir das ungeordnete XY-Modell im Bereich 

rj-i ^ rp* 



Fiir die Fugazitat y in (4.39) ist der zum jeweiligen Gittertyp passende Ausdruck 
anzusetzen. Offenbar entspricht Gleichung (4.39) der in Abschnitt 2.3.2 gewonncncn 
Gleichung (2.59). Dort konnte bereits gezeigt werden, dafi der hieraus mittels 
^ = ~WM folgende Flufi der Entropie im Fall T > T* positiv ist. Beim 
Ubertragen der dort benutzten Argumente auf ( 4.39[ ) ist lediglich a durch j zu 
ersetzen. Fiir T > T* liegt somit kein unphysikalisches Verhalten vor. 

Der interessantere Fall ist T < T* . Hier wurde der Entropieflufi bei dem von 
Nelson vorgeschlagenen Verhalten negativ und hier kommen auch die in dieser 
Arbeit entwickelten Modifikationen zum Tragen. Fiir das Quadratgitter folgt in 
diesem Temperaturbereich aus (4.37) unter Benutzung von (3.50) 



d(dFy) 



= —irTrdr 



1 



e -PEr 



2ti 



dO 
— < 

2tt 



(4.40) 
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wobei wie in Anhang B hergeleitet A 2 = ^£-(hir — cos 2 6) gilt. Integration 
liefert als wesentlichen Term auf grofien Langenskalen 

f 27T d6 A 2 (£y+g m ) 2 
dF v = -nrdr / e ^ . (4.41) 

Setzt man hier die expliziten Ausdriicke fur £V, E w und A ein, so findet man 
auf geniigend grofien Langenskalen (lnr 3> ^f- und \nr ^> a) 

dF v = -T^nT* rdrP{r), (4.42) 



und fur den Flufi der freien Energiedichte im Bereich T < T 

-T*r 2 P{r) 



dFy n m* 2 



dl 

(4.43) 

Auch dieses Ergebnis stimmt mit dcm von Tang und Schcidl gcfundenen Verhalten 
des ungeordneten XY-Modells bis auf Vorfaktoren iiberein. Fiir das Dreiecksgitter 
ergibt sich bei Verwendung der passenden Wahrscheinlichkeit P(r) und der 
passenden Fugazitat y die gleichc Beziehung. 

Zur Be rechn ung des Entropieffusses mittels ^f- — —-§f^[f- beschranken wir 
uns in ( |4.43| ) auf die wesentlichen Faktoren, die das Verhalten auf grofien Skalen 
bestimmen. Wir betrachten daher nur ^jj- oc ~T*y 2 . Die vernachlassigten 
Faktoren haben keinen Einflufi auf den Entropicflufi. 

Wegen 

|(n 2 ) = 0, (4.44) 

was sowohl beim Quadrat-, als auch beim Dreiecksgitter gilt, liegt kein Entropieflufi 
vor. Zwar wurden hicrbei die impliziten Temperaturabhangigkeiten vernachlassigt, 
da allerdings auch 

^4(7V)=0 + O(y 4 ) (4.45) 

fiir beide Gittertypen gefunden wird, liefern auch diese impliziten Abhangigkeiten 
(vgl. hierzu auch Abschnitt 2.3.2 und j(| ) keinen Entropieflufi fiir T < T*. 
Dieses Verhalten, welches Tang || ebenso fiir den Tieftemperaturbereich des 
ungeordneten XY-Modells findet, erklart er mit einem Einfrieren der Defektpaare 
bei T = T* . Festzuhalten bleibt, dafi anders als bei Nelson || der Entropieflufi 
im Bereich T < T* nicht unphysikalisch wird. 
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Kapitel 5 

Einflufi eines Substrats 



Bishcr habcn wir lediglich cincn frcicn zwcidimcnsionalcn Film bctrachtet. Unter 
realen Bcdingungcn licgt cin solcher Film allerdings immer auf einem Substrat 
auf, mit dem dcr Film auch wechsclwirkt. Die in den vorangegangenen Kapiteln 
gefundenen Ergebnisse sind in diesem Fall nur dann anwendbar, wenn man 
die Wechselwirkung mit dem Substrat auf groBen Langenskalen vernachlassigen 
kann, die Ordnung im kristallinen Film also nicht gestort wird. 

Fiir ein glattes Substrat ist dies erfullt, da die Wechselwirkung zwischen Substrat 
und kristallincm Film konstant ist. In Expcrimcnten ist ein Film aber meist 
auf Substrate aufgebracht, die, wie der Film selbst, eine kristalline Struktur 
aufweisen, und deren Gitterkonstante die gleiche Grofienordnung besitzt, wie 
die Gitterkonstante des Films. 



5.1 Kommensurables Substrat 

Wir beschranken uns hier auf die Untersuchung einer kommensurablen Situation, 
bei der das Substrat ein Potential erzeugt, dessen Minima dieselbe Geometrie 
wie der kristalline Film aufweisen, das heiBt im Fall eines kristallinen Films mit 
Quadratgitter-Struktur sind die Minima ebenso als Quadratgitter angeordnet, 
und im Fall eines Dreiecksgitters bilden auch die Potentialminima cin Drciecksgitter. 
Ist a die Gitterkonstante des kristallinen Films, so gelte fiir die Gitterkonstante 
b des aus den Potentialminima gebildeten Gitters 

&=-, (5-1) 
n 

wobei n eine natiirlichc Zahl sei. Liegt ein solches Potential vor, ist es moglich 
allc Gittcrbaustcine des kristallinen Films in einem Potentialminimum zu plazieren, 
ferner ist kein Potentialminimum mehrfach besetzt. Es ist wesentlich zwischen 
der Gitterkonstanten b des aus den Potentialminima gebildeten Gitters und der 
Gitterkonstanten as des Substrats zu unterscheiden, da b ^ as gelten kann. Ein 
Bcispicl hierfiir zeigt Abbildung 6.2, wo die Potentialminima ein Dreiecksgitter 
bilden, die das Potential erzeugenden Substratatome konnen allerdings bei den 



•57 



5.1. KOMMENSURABLES SUBSTRAT 



Potentialmaxima sitzen. In diesem Fall waren sie in einem hexagonalen Gitter 
mit Gitterkonstanten = angeordnet. 

Fiir T = und a = werden daher alle Kristallbausteine in den Potentialminima 
eingelagert sein und ein maximaler Energicgcwinn vorliegen. Erst bei hoheren 
Temperaturen oder starkerer Unordnung wird es moglich sein, die Wechselwirkung 
zwischen Substrat und Film auf groBen Langenskalen zu vernachlassigen. 

Im Fall des Quadratgitters sind die primitiven Gittervektoren des kristallinen 
Films 

*-(;) *-(;), (5. 2) 

fiir die primitiven Vektoren des reziproken Gitters folgt hieraus 




Diese Vektoren erfiillen die Relation bi Gj = 2nSij . Analog hat man fiir das 
Dreiecksgitter 




und 

Gx = ( i_ ) G 2 = ( 4° V (5.5) 



Nun soil die periodische Wechselwirkungsenergie V^(u) eines Kristallbausteins 
mit dem Substrat in Abhangigkeit der Verschiebung u aus seiner Gleichgewichtsposition 
in einem Potentialminimum dargestellt werden. Eine geeignete allgemeine Formulicrung, 
die die oben aufgefiihrten Bedingungen an das Potential erfullt, liefert die Darstellung 
in Form einer Fouricrrcihc 

V(u) = ~Y^ Cfcj exp(i n(kGx + lG 2 )u). (5.6) 

fc,; 

Hierbei sind k und I ganze Zahlen. Fiir den Vorfaktor mufi Cm — C-k-i 
gelten, damit das Potential V reell bleibt, ferner fordern wir Cki > 0, um 
Potentialminima garantieren zu konnen. Setzt man Gki = n(fcGi + ZG2), so 
ist Gki immer ein Vektor des reziproken Gitters. Es ist also 

V(u) = -J2 C ^ ex P( l G ki u )> ( 5 - 7 ) 

fc,; 



die Summation entspricht hierbei einer Summation iiber alle reziproken Gittervektoren. 
Zur beispielhaften Darstellung solcher Potentiale beschranken wir uns auf jene 
Summanden aus ( |5.7| ), bei denen Gki einen minimalen, nicht verschwindenden 
Betrag besitzt. Wie (5.17) zeigen wird, tragen gerade diese Terme wesentlich 
zur Hamiltonfunktion bei, Terme mit grofieren Betragen von Gki konnen fiir 
die Fragestellung dieses Kapitels vernachlassigt werden. 
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Fur das Quadratgitter haben (±Gi = ±nGi) und (±Goi = ±nG 2 ) minimalcn 
Betrag. Setzen wir mit C\q = Coi = v ein beziiglich der Rotation inn ^ 
symmetrisches Potential an, so resultiert 



V(u) = — v [cos (n Giu) + cos (n G2U)] 



(5. 



2tt 

cos I n — u a 
a 



2ir 

cos I n — u 

a 



(5.9) 



Dieses Potential wird in Abbildung 5.1 gezeigt. 

A separate downloading of this figure is possible. 

Abbildung 5.1: Einfaches Potential einer Wechselwirkung zwischen einem Kristallbaustein 
und dem Substrat fiir das Quadratgitter. Dunkle Regionen entsprechen Potentialminima, 
helle den Potentialmaxima 

Fiir das Dreiecksgitter lauten die Gki mit minimalem Betrag (±Gio = ±nGi), 
(±G i = ±nG2) und (±G_n = ±n(G 2 — Gi)). Die analoge Darstellung des 
Potentials (vgl. Abbildung 5.2) ist: 



V(u) = — v [cos (n Giu) + cos (n G 2 u) + cos (n (G2 — Gi)u) 



(5.10) 



47T 

n — = u n 



a 

2tt 
-n — u x 
a 



2tt 
n — u x 
a 



2ir 

■ n — 7= u, 



iV5 



2tt 
a\/3 



(5.11) 



A separate downloading of this figure is possible. 

Abbildung 5.2: Einfaches Potential einer Wechselwirkung zwischen einem Kristallbaustein 
und dem Substrat fiir das Dreiecksgitter. Dunkle Regionen entsprechen Potentialminima, 
helle den Potentialmaxima 



Im folgenden kehren wir wieder zur allgemeinen Formel (5.7) zuriick. Betrachtet 
man nun nicht nur die Wechselwirkung eines Kristallbausteins mit dem Substrat, 
sondern die Wechselwirkung aller Bausteine, so ergibt sich in der Kontinuumsform 
folgender Beitrag zur Hamiltonfunktion 



Hs U b = -v C, 

kl 



kl 



d 2 r exp(zG k iu(r)). 



(5.12) 



Hierbei ist v := max{Cfc;}. Die Cm aus (5.7) miissen ein solches Maximum 
haben, da sonst ein Kristallbaustein durch die Wechselwirkung mit dem Substrat 
beliebig viel Energie gewinnen konnte, was unphysikalisch ist. Ferner gilt Cm '= 
, woraus < Cm < 1 folgt. 

Die Kontinuumsform in ( |3.12 ) ist gerechtfertigt, wenn ein kontinuierliches Verschiebungsfeld 
u(r) vorliegen kann. Dies ist bei sehr groBen Werten von v nicht mehr der Fall. 
Die starke Bindung der Kristallelemente in den Potentialminima wurde dazu 
fiihren, dafi u(r) nur noch diskrete Werte annehmen kann, da jede Kristallelement- 
Verschiebung u ein Gittervektor des Potentialminima-Gitters sein miifite. 
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Einc Bcdingung fur v lafit sich wiefolgt abschatzen. Fiigt man einem Kristall, der 
kommcnsurabel (mit n = 1) auf einem Substrat aufliegt, eine zusatzliche Reihe 
von Kristallbausteinen hinzu, so lagert sich diese Reihe im Falle sehr grofier 
Werte von v bei einer anderen, bereits in Potentialminima liegenden Reihe an. 
Diese Minima waren somit doppelt besetzt, wahrend links und rechts davon die 
normale kommensurable Situation erhalten bliebe. Der Storbereich, den man 
auch Domanenwand nennt, hatte somit eine Breite von Iq = b. Im Falle kleinerer 
v hingegen wurde sich die Breite der Domanenwand vergrofiern. Ein Vergleich 
der elastischen Energiekosten aufgrund von u(r) und des Energiegewinns durch 
die Wechselwirkung mit dem Substrat zeigt, dafi ein Energicminimum bei einer 
Domanenwandbreite von 



i^ b .m^±. (5.i3) 



vorliegt . Da die Bedingung fur die Kontinuumsform Iq ^> b lauten mufi, 
folgt fur die Starke des Potentials v -C ^rp^ ■ 

Der Beitrag von Hs u b lafit sich dann als Storung der Hamiltonfunktion auffassen. 
Storungsrechnung in erster Ordnung liefert 

Hsub = ~vJ2 Ckl f d 2 r [<exp(iGidu(r)))] D , (5.14) 

kl J 

wobei < ... > eine thermische Mittehmg mit dem ungestorten Hamiltonian und 
[■■■]d eine Mittclung iiber die Unordnung symbolisieren. Da in der kristallinen 
Phase auf sehr grofien Langenskalen die Versetzungen keine Rolle mehr spielen, 
ist eine Mittclung mit einem rein phononischen Verschiebungsfeld moglich, welches 
eine Gaufische Wahrscheinlichkeitsvcrtcilung zeigt. Berechnungcn von Nelson § 
verwendend (vgl. Anhang D), finden wir 

_ TG 2 (3m + A) _ gG 2 njj(n + A) 2 

[(exp(iGu(r)))] = R R ^^+>-r 2 . (5.15) 

Hierbei sind R die Systemgrofie, Slo die in Kapitel 2 erlauterte Anderung der 
Kristallnlmnache aufgrund von Verunreinigungen und A, fi die vollstiindig renormierten 
Lame-Koefhzienten, die zu verwenden sind, da obige Rechnung ja sehr grofie 



Langenskalen bedingt. Betrachtet man ( 5.14 ), so stellt man fest, dafi der Beitrag 
des Summanden kl mit 

8vr M (A + 2 M ) 8^(2^ + A) 2 1 ' ' 

skaliert (d ist hier die Dimension). Der elastische Anteil der Hamiltonfunktion 
skaliert wegen der dort enthaltenen Ableitungen mit d — 2. Im vor liegenden 
zweidimesionalen Fall (d = 2) folgt daraus, dafi bei grofien Systemen der Summand 
kl somit nur dann wesentlich zur Hamiltonfunktion beitragt, wenn 

TG 2 fci (3 M + A) aGj^ + Xf 
~ 8^(A + 2^) 8^(2^ + A) 2 ' 1 ' ' 

erfiillt ist. 
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1st diese Ungleichung fur die Summandcn mit den kleinsten Betrage von Gki 
verletzt, so gilt sic fur kcincn Summanden. Daher sind nur diese Summandcn 
fur unsere Betrachtung wesentlich (vgl. Bemerkung vorne). Beim Quadratgitter 
liefern die Gki mit den kleinsten Betragen (Gio und Goi) G 2 = Gqj = n 2 ^-, 
fur das Dreieckgitter gilt analog G\ = Gq X = G?_ u = n 2 ^^-. 



Eine Verletzung der Ungleichung ( |5.17| ) fiir alle Summanden bedeutet, daB 
die Wechselwirkung mit dem Substrat vernachlassigt werden kann. Die in den 
vorangehenden Kapiteln untersuchte Theorie des Schmelzens ist dann weiterhin 
anwendbar. Die Bedingung fiir die Existenz eines Kosterlitz-Thouless-Phasenubergangs 
vom kristallinen Zustand mit quasi-langreichweitiger Ordnung in den nicht- 
kristallinen Zustand lautet im Falle eines unterliegenden kommensurablen Substrats 
daher: 

T(3 M + A) anKfi + X) 2 j jjgp fiir Quadratgitter 

\ (in 2 ■■■■ ■ , ■.. K^ 1 *) 



8tt/i(A + 2fi) 8tt(2^ + A) 2 1 T J^- T fiir Dreiecksgitter 

Dieses Ergebnis wurde fiir Systeme ohne Unordnung bereits von Pokrovsky und 
Talapov j29| gefunden. 

Kritischer Wert von n 

Im folgenden betrachten wir lediglich das Dreiecksgitter, da die benutzte Elastizitatstheorie 
exakt nur fiir diese Gitterform gelten und man die analogen Ergebnisse fiir das 
Quadratgitter durch Austausch von konstanten Faktoren erhalt. 



Es ist nun notwendi g die Bedingung ( 5.18 ) mit der in Abbildung 4.5 dargestellten 



und durch Formel (4.12) beschriebenen Phasengrenze zu vergleichen. Mit der 



Kopplungskonstanten 

J= 4 -t 

2/J. + X 



j= 4 ^ i+ y ( 5.i9) 



und unter Verwendung von a — a— f wird (5.18) zu 



T (3 M + A)(/j + A) _ Ox + A) 2 3 

J (2 M + A) 2 4(2/i + A) 2 47rn 2 ' [ 1 



Beachtet man die physikalischcn Bedingung, daB der Schubmodul fi und der 
Kompressionsmodul \x + A beide positiv sein miissen [^3| , so findet man, daB die 

Grofien (^^jt^ un d ^+x)' 2 aus ( 5.20| ) nur Werte inner halb des Intervalls 



[0, 1] annehmen konnen. Das Maximum 1 wird im Fall A 3> (J. erreicht. Da ein 
unterer kritischer Wert von n berechnet werden soil, sind diese beiden Werte 
maximal zu setzen. Somit gilt 

H>iJ?- < 5 - 2i > 

Der kritische Wert von -j, bis zu dem die kristalline Phase noch existieren 
kann ist der entsprechende kritische Wert der Unordnungsstarke ist a c — 
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yg^r. Um die kristalline Phase mit QLRO auch bei unterliegendem Substrat 
beobachctcn zu konncn, mufi 



1 

16tt 



1 



> 



64-7T 47m 2 

erfullt sein, was der notwendigcn Bcdingung n > 4 entspricht. 



(5.22) 



Es bleibt noch zu zeigen, dafi n > 4 auch eine hinreichende Bedingung ist, das 
heifit, da fi es im Fall A ^> /i tatsachlich einen Bereich der geordneten Phase gibt 
fur den ( |5.1S| ) fur 
moglich) . 



4 erfullt ist 



4- und er 

1D7T 



16 



g^r sind ja nie gleichzeitig 



Nur im System ohne Unordnung liegt der Phaseniibergang beim Maximalwert 



Yg^-. (5.21) reduziert sich hier zu 



T 3 
J 47m 2 ' 



(5.23) 



Fiir n = 4 ist diese Ungleichung erfullt, somit ist n > 4 auch eine hinreichende 
Bedingung fiir die Existenz der kristallinen Phase mit QLRO. 

Abbildung 5.3 verdeutlicht dieses Ergebnis, indem die Bedingung ( 5.20| ) graphisch 
gezeigt wird. Hierbei wird ein System ohne Unordnung mit X = 2fi angesetzt. 
Man sieht deutlich, dafi nur fiir n > 4 die kristalline Phase mit QRLO existiert. 



kristallin 


ungeordnet 


kommen- 




surabel 






T 



Abbildung 5.1: Beispielhafte graphische Darstellung von (5.2C ) fiir \ = 2fi und a — 



Insgesamt ist also festzustellen, dafi bei einem unterliegenden kommensurablen 
Substrat, der beschriebene Kosterlitz-Thoulcss-Schmclzubergangnur dann gefunden 
werden kann, wenn die Gitterkonstante des Substrats die Relation b < f erfullt. 
Ist dies nicht der Fall, fallt die gesamte kristalline Phase mit QLRO weg, da sie in 
dem Bereich lage, in dem die Wechselwirkung mit dem Substrat auch auf grofien 
Langenskalen wichtig ist. Daher tritt auch der Kosterlitz-Thoulcss-Ubergang 
von der kristallinen in die hexatische Phase nicht mehr auf. Als kritischen Wert 
fur n haben wir somit 

n knt = 4 (5.24) 
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Fiihrt man obige Betrachttmgen vollig analog fur das Quadratgitter durch, 
so findet man denselben kritschen Wert fur n. Das Ergebnis deckt sich mit 
der qualitativen Aussage von Lyuksyutov, Naumovets und Pokrovsky, daB die 
kristalline Phase mit quasi-langreichweitiger Ordnung nur fur grofiere Werte von 
n gefunden werden kann jl6| . Ziel des folgenden Abschnitts ist es, bei festem 
n ein Phasendiagramm fur einen kristallinen Film zu finden, der auf einem 
kommensurablen Substrat aufgebracht ist. 



5.2 Phasendiagramm bei unterliegendem Substrat 

Sowohl der Schub- als auch der Kompressionsmodul sind positiv. Dies entspricht 
fur die Lamc-Koeffizicntcn den Bcdingungcn /i > und —\i < A < oo. Setzen 
wir 

k = - (5.25) 
A 1 

so liegt die Konstante k im Intervall [—1, oo]. Messungen an rcalcn physikalischen 
Systemen zeigen, dafi die Werte fur n tatsachlich in einem Intervall [—1, 10] zu 
finden sind. 

Mochte man bei festem n den Bereich der geordneten Phase des in Kapitel 
4 gewonnen Phasendiagramms abschatzen, der bei unterliegendem Substrat 
noch beobachtbar ist, so ist es sinnvoll diese Abschatzung in Abhangigkeit der 



Konstanten k durchzufiihren. Fiihrt man n in (5.20) fiir das Dreiecksgitter ein, 
so erhalt man: 

T (3 + «)(! + «) _ {l + nf 3 

J (2 + K ) 2 + 4(2 + K ) 2 W [ ' 

Fiigt man dem Phasendiagramm aus Abbildung 4.5 noch eine dritte K-Achse 

hinzu, so lassen sich sowohl die Phasengrenze des Kosterlitz-Thouless-Schmelziibergangs, 



die keine K-Abhangigkeit zeigt, als auch die Bedingung (5.26) als Flachen in 
diesem Phasendiagramm darstellen. So kann der Bereich des Phasendiagramms 
verdcutlicht werden, in dem die Wechselwirkung mit dem Substrat wichtig wird 
und die in Kapitel 4 erhaltenen Resultate daher nicht mehr zutreffen. 

Im folgenden wird das Phasendiagramm fiir das Dreiecksgitter bei n — 6 dargestellt 
und erlautert. 



Abbildung 5.4 zeigt die Abhangigkeit der kritischen Temperatur von k bei 
(7 = 0. Da, wie bereits erlautert, der Wert fiir ^ des Kosterlitz-Thouless- 
Phasenubergangs keine k- Abhangigkeit besitzt, liegt die Phasengrenze zwischen 
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kommen- 
surabel 



kristallin 



ungeordnet 



Abbildung 5.2: Phasendiagramm 
von k ohne Unordnung (a = 0). 



Dreiecksgitters bei n — 6: -j in Abhdngtgkeit 



fliissiger Phase ohne Translationsordnung und fester kristalliner Phase mit quasi- 
langreichweitiger Translationsordnung konstant bei j^. Die andere eingezeichnete 
Kurve separiert den Bereich, in dem (5.18) erfiillt ist, von jenem in dem diese 
Bedingung nicht gilt. Rechts dieser Kurve (hier gilt ( 5.18 )) spielt die Wechselwirkung 
mit dem Substrat auf grofien Langenskalen somit keine Rolle. Hier existieren die 
kristalline Phase mit QLRO und der Kosterlitz-Thouless-Ubergangzur ungeordncten 
Phase. Im linken Bereich ist die Wechselwirkung mit dem Potential wesentlich. 
Wir erwarten hier, dafi die Gitterbausteine in die Potcntialminima eingebettet 
sind und langreichweitige Ordnung zeigen. Es liegt also eine kommensurable 
Phase vor. Somit bildet die Bedingung (5.18) die Phasengrenze zwischen der 
kommensurablen Phase mit langreichweitiger Ordnung und der ublichen kristallinen 
Phase mit quasi-langreichweitiger Ordnung. Dieser Phasenubergang ist kontinuierlich. 
Da (5.15) in der ungeordneten Phase nicht gilt, kann nicht erwartet werden, dafi 
die Phasengrenze der kommensurablen Phase in diesen Bereich fortsetzbar ist. 
Sie ist daher dort nur gestrichelt eingezeichnet. Im schraffierten Bereich ist ein 
Phasenubergang von der kommensurablen zur ungeordneten Phase zu erwarten, 
dessen genaue Lage und Ordnung hier nicht diskutiert werden soil, ebenso wie 
die Frage, ob in diesem Bereich eine mogliche Zwischenphase existiert. 

Analog zur vorhergehenden Abbildung zeigt Abbildung 5.5 nun den Phasenverlauf 
in einem cr-K-Diagramm. Hier ist T = gj- fest gewahlt. Man sieht, dafi auch hier 
bei geniigend niedriger Unordnungsstarke eine kommensurable Phase existiert, 
die bei niedrigen Werten von k in die ungeordnete Phase iibergehrt, wahrend 
bei steigendem k die kristalline Zwischenphase zu finden ist. 

Darstellung 5.6 zeigt das er-T-Diagramm fur ein festes Verhaltnis zwischen den 
Lame-Koeffizienten k — 2, welches von realen Systemen durchaus angenommcn 
werden kann. Dieses Diagramm ist mit Abbildung 4.5 zu vergleichen. Die kommensurable 
Phase, die bei tiefen Temperaturen und geringer Unordnungsstarke auftritt, ist 
durch eine gerade Phasengrenze von der kristallinen Phase separiert. Bei sehr 
niedrigen Temperaturen geht sie direkt in die ungeordnete Phase iiber. Lchncn 
|30[ hat ein solches Phasendiagramm fur ein ungeordnetes XY-Spinsystem 
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Abbildung 5.3: Phasendiagramm des Dreiecksgitters bei n = 6: a in Abhdngigkeit von 
n bei gegebener Temperatur T — . 



gefunden, welches nicht durch ein Substrat, sondern durch ein aufieres Magnetfcld 
gestort wird. Neben den auch hicr gefundenen Phasen mit vergleichbarcn Phascngrenzen, 
diskuticrt cr im schrafficrtcn Bcrcich die Existenz einer zusatzlichen superrauhen 
Zwischenphase (In 2 r-Abhangigkeit der Korrelationsfunktion), deren Grenzen 
aber nicht genau bestimmt werden konncn. 

Festzuhalten bleibt, da8 sich die Grenze der kommensurablen Phase mit steigendem 
n nach links verschiebt, der schrafherte Bereich entfallt dann. Bei nicdrigcrcn 
Werten von n allerdings verschiebt sich diese Grenze nach rechts. Fur n < 4 
befindet sie sich soweit rechts, dafi die kristallinc Phase nicht mehr auftritt. Hier 
cxisticrt dann der vornc beschricbenc Kostcrlitz-Thouless-Schniclziibergang nicht 
mehr, sondern lediglich die kommcnsurable und die ungcordncte Phase. Somit 
ist insbesondere fur a = b kein Kosterlitz-Thouless-Schmelzubergangzu erwarten. 

Abbildung 5.7 zeigt abschliefiend nocheinmal eine dreidimensionale Darstellung 
der oben diskutierten Phascngrenzen. Im linkcn Bereich bei klcincn Werten von 
d und -y befindet sich die kommensurable Phase, wahrend im rechten Bereich, 
bei hohen Werten von ^ und k, die ungeordnete Phase vorliegt. Die kristalline 
Phase mit quasi-langreichweitiger Ordnung befindet sich zwischen den beidcn 
Phascngrcnzflachen. 

A separate downloading of this figure is possible. 

Abbildung 5.7: Dreidimensionales Phasendiagramm des Dreiecksgitters bei n = 
6: Hier sind die Unordnungsstarke a, -j und n gegeneinander aufgetragen. Die 
Fortsetzung der Phasengrenze zwischen kommensurabler und kristalliner Phase 
in den Bereich der ungeordneten Phase hinein ist, wie bereits erwahnt, nicht 
korrekt, aus Griinden der Anschaulichkeit hier aber dennoch eingezeichnet. 

Neben den beschriebenen kommensurablen Situationen, in denen die primitiven 
Gittervektoren des aufliegenden Films und des Potentialminimagitters bis auf 
einen Faktor n gleich sind, existieren auch kommensurable Situationen, bei 
denen dies nicht der Fall ist. Abbildung 5.8 zeigt ein solches Beispiel. Hier ist 
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Abbildung 5.4: Phasendiagramm des Dreiecksgitters bei n — 6: Unordnungsstarke a 
in Abhangigkett der Temperatur, k = 2 ist hier fest gewahlt. 



das Drciccksgitter der Potcntialminima gegeniiber dem Dreiecksgitter des Films 
(weiBe Punkte) um 30 Grad gedreht. 

A separate downloading of this figure is possible. 

Abbildung 5.8: Potential der Wechselwirkung zwischen einem Film (weifie Punkte) 
und einem kommensurablen Substrat. Die beiden Dreiecksgitter sind um 30 Grad 
zueinenader gedreht. 



Hier ist a = V3b. Mochte man die Situation analog zu (5.1) verallgemeinern, so 
gilt 

a = nV3b. (5.27) 



Geht man wie in Abschnitt 5.1 vor, so findet man anstelle der Bedingung (5.18) 
in diesem Fall: 

T(3^ + A) q^ + A) 2 2a 2 

> 1G -2„2 - I 5 ' 28 ) 



8tt/z(A + 2/j.) 8tt(2^ + A) 2 
Als kritischer Wert ergibt sich daraus nkrit — 2. 

In analoger Weise lassen sich alle kommensurablen Situationen behandeln, bei 
denen sich nicht mehr als ein Kristallbaustein in einem Potentialminimum befindet. 

Nicht diskutiert wurden in diesem Kapitel 

• Kommensurable Substrate hoherer Ordnung mit einer Gitterkonstanten 

b = pa, (5.29) 

wobei p eine natiirliche Zahl ist. Hier befinden sich dann p Kristallbausteine 
in einer Potentialmulde. Nach |2jJ wird in diesem Fall die rechte Seite 
der Ungleichung ( |5.18 ) proportional zu p-. Insgesamt ist also wiederum 



eine Proportionaltat zu a 2 festzustellen, wie schon im Fall b = a. Es kann 
also bei diesen Substraten ebenfalls kein Kosterlitz-Thouless-Schmelzen 
auftreten. 
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• Inkommensurable Substrate mit einer Gitterkonstanten der Form 

L P , 



n 



■{a -6), (5.30) 



hierbei ist 8 = der Misfit zwischen den beiden Gitterkonstanten. 
Entsprechend |n] ist mit steigendem Misfit ein Absenken der Phasengrenze 
der kommensurablen Phase zu niedrigeren Temperaturen und Unordnungsstarken 
zu erwarten, bis zu einem kritischen Wert von 8, bei dem die kommensurable 
Phase vollstandig verschwindet. 

• Substrate mit beliebiger Gitterkonstante, die aber gegeniiber dem Gitter 
des kristallinen Films gedreht sind; die Symmetrieachsen der beiden Gitter 
sind also nicht parallel. Sofern eine kristalline Phase mit quasi-langreichweitiger 
Ordnung und damit ein Kosterlitz-Thoulcss-Schmelzubergang in die hcxatischc 
Phase existieren, findet man in ( 2.42] ) zwei verschiedene Kopplungskonstanten 



fur den Logarithmus-Term und den Winkelterm, die unterschiedlich renormieren 
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Kapitel 6 

Zusammenfassung 



In dieser Arbeit wurde das Schmclzcn zwcidimensionaler kristalliner Filme bei 
vorlicgendcn eingefrorenen Verunrcinigungen untersucht. Im Rahmen der Koster- 
litz-Thouless-Theorie wird ein zweistufiges Schmelzen mittels zwei kontinuierlichcr 
Phaseniibergange erwartet. Zunachst geht die kristalline Phase mit quasi-langreichweitiger 
Translationsordnung (algebraische abfallende Translations-Korrelationsfunktion) 
in eine hexatische bzw. tetratische Flussigkristall-Phase iiber, die keine Translationsordnung 
mehr aufweist, aber noch quasi-langreichweitige Orientierungsordnung zeigt. 
Beim zweiten Phaseniibergang in die isotrop fLiissige Phase geht auch diese 
verloren. Gegenstand dieser Arbeit ist lediglich der erste Phaseniibergang in 
die Flussigkristall-Phase. 

Dieser ist bereits von Nelson [^| im Falle glatter Substrate betrachtet worden. Er 

fand bei tiefen Temperaturen und Unordnungsstarken a > einen Wiedereintritt 

in die Flussigkristall-Phase ohne Translationsordnung. Im Bereich T < 4^, der 

den gesamtcn Wiedereintritt umfafit, fiihren die von ihm gefundenen RenormicrungsfluBglcichungen 

der Kopplungskonstanten J und der Fugazitat y allerdings zu einem Absinkcn 

der Entropie bei Renormierung und zu negativen Entropiewerten auf groBen 

Langenskalen. Grund hierfur war, dafi Nelson fur die Wahrscheinlichkcit einen 

Versetzungsdipol vorzufinden eine Boltzmann-Verteilung angesetzt hat, die im 

Bereich des Wiedereintritts sehr grofie Werte annimmt. Damit wird dort aber 

auch die Fugazitat sehr groB und die iibliche Naherung fiir kleine Fugazitaten 

kann nicht mehr angewendet werden. 

In dieser Arbeit wird eine normierte Dipolwahrscheinlichkeit angesetzt, die 
Versetzungen als nicht wechselwirkende Fermionen behandelt. Dies geschieht 
in Anlehnung an eine Arbeit von Nattermann, Scheidl, Li und Korshunov Q 
, die dieses Konzept fiir das ungeordnete XY-Spinsystem entwickelt haben, 
wo Vortices als topologische Defekte auftreten. Renormierung mit Hilfe eines 
dielektrischen Formalismus liefert FluBgleichungen fiir die Kopplungskonstante 
und die Fugazitat, die fiir T > 4p mit den von Nelson gefundenen Gleichungcn 
ubereinstimmen. Eine Renormierung der Unordnungsstarke wird nicht untersucht. 
Das aus den Flufigleichungen resultierende Phasendiagramm zeigt keinen Wiedereintritt, 
im Bereich T < 4p verlauft die Phasengrenze zwischen fester Phase und Flussigkristall- 
Phase parallel zur Temperaturachse bei einem Wert von a c = j^-. Abbildung 
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4.5 zeigt das Phasendiagramm noch einmal. Dieses Resultat stimmt mit den 
Abschatzungen der Arbeit von Cha und Fertig M iiberein. Ein Absinken der 
Entropie bei Renormierung wird nicht mehr gefundcn. 

Zusatzlich wurde der Einflufi eines periodischen Substratpotentials nicht-glatter 
Substrate untersucht. Diese Analyse beschrankte sich auf kommensurablc Situationen, 
bei denen die Potentialminima den selben Gittertyp bilden wie der aufliegende 
kristallinc Film. Hat der Film die Gitterkonstante o und das Gitter der Potentialminima 
die Gitterkonstante 6, so wurde a = nb gcfordert, wobei n eine natiirlichc 
Zahl ist. Es konnte gezeigt werden, dafi die Wechselwirkung mit dem Substrat 
nur dann die Ausbildung einer kristallinen Phase mit quasi-langreichweitiger 
Translationsordnung zulafit, wenn n > 4 erfiillt ist, das Substrat somit schnell 
moduliert (vgl. Abbildung 5.3). Anderenfalls, also auch fur a = b, verschwindet 
diese Phase vollig und damit auch der oben beschriebene Phaseniibergang. An 
ihre Stelle tritt eine kommensurable Phase mit langreichweitiger Ordnung, die 
direkt schmilzt. Fur n > 4 existiert diese kommensurable Phase zwar auch, 
geht aber zunachst in die kristalline Phase mit quasi-langreichweitiger Ordnung 
liber, bevor dann der in dieser Arbeit behandelte Phaseniibergang eintritt (vgl. 
Abbildung 5.5). 

Ein Vergleich mit Experimenten und Simulationen bringt keine direkte Bestatigung 
der Ergebnisse dieser Arbeit, da sich sowohl Experimente als auch Simulationen 
bishcr auf die Frage nach der Richtigkeit der Kosterlitz-Thouless-Beschreibung 
beschranken und daher noch keine Arbeiten vorliegen, die versuchen den Einflufi 
von Unordnung oder verschiedenen Substratpotentiale auf den Schmclzvorgang 
explizit zu untersuchen. Dennoch scheint das Verhalten von adsorbiertem Krypton 
auf einem kommensurablen Graphitgitter die hier erzielten Resultate zu bestatigen. 

Mit dieser Arbeit konnte somit das in jiingster Zeit gefundene Verhalten des 
ungeordneten XY-Modells auf das verwandte Modell zweidimensionaler Gitter 
mit eingefrorenen Verunreinigungen ubertragen werden. 
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Anhang A 

Coulombgas-Beschreibung 



A.l Abbildung auf ein Coulombgas 



Ausgehend von ( |2.4C| ) und unter Verwendung der Formeln (2.36) und (2.39) 
erhalt man: 

H e i = / / d 2 rd 2 r'h{T)e l0 -,V M^kN' l (r - r') 2 (In |r - r'| + C) 



16tt 



|r-r'|>l 

(A.l) 



Benutzt man b(r) = ^ Q b a (5(r — r a ), so resultiert mit einer modifizierten von 
Nelson [e3 benutzten Konstanten C — C + \ 

H e i = -j^Y^ajtpbaieij V a jbj3 k €ki(x a i - x i) (hi\r a - ra| + C\ 

= b ai b 0i (in KfJ \ +C) + WEHi^glW^^ 



J 



E«# b a ib/3i (in |r a/3 | + Cj + & ai 6/3i 



(f)e»i(^ a i-Z0i))(fc0j(g e «j-3:gj)) 

Unter Benutzung der Ladungsneutralitat ^h a = gilt die Beziehung 

-(<7 + l)^b a b /3 = (C' + l)^b a 2 . (A.3) 

Dcfiniert man nun Ec — (C + l)J/8w als Energie eines Versetzungskerns, so 
erhalt man direkt Formel ( 2.42|) , welche zu zeigen war: 



8?r 



H el = -4zY. fbab^lnlr^l - + ^ |b Q | 2 (A.4) 



Dies ist exakt das Ergebnis von Nelson [g, |24| ■ Young [jl3| hat eine leicht 
andere Definition von Ec benutzt und daher mit einem zusatzlichen Term in 
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dcr Hamiltonfunktion gearbeitet. Wir folgcn im weiteren aber dem Formalismus 
von Nelson. C ist eine positive Konstante, deren Wert von der Gitterstruktur 
abhangt, sic mi8t das Verhaltnis von Kerndurchmesser zu Gitterkonstante. Ihr 
Wert ist von der Ordnung 0(1), fur das Quadratgitter wurde C w ^ berechnet. 
Da fur die Erzeugung eines Versetzungsdipols mit Separation 1 (a = 1: Gitterkonstante) 
maximal die Energie 2E C aufgewendet werden mufi, kann E c tatsachlich als 
Energie eines Versetzungskerns bezeichnet werden. 
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A. 2 Vergleich mit der Elektrodynamik 

In diesem Abschnitt sollen die in dicscr Arbeit auftretenden 
elastizitatstheoretischen Grofien mit den analogen elektrodynamischen Groficn 
verglichen werden. Die folgende Tabelle zeigt die wesentlichen Relationen 



Elastizitatstheorie 

Quellenfunktion 



r)(r) = tijVi ^2 b a,% S(r - r„) 



Elektrodynamik 

Ladungsdichte 



e( r ) = 53 Sa6( - r ~ 



Greenfunktion 



1 



g(r) = —r In |r| 



Greenfunktion 

g(r) = ln|r| 



Verzerrungsfunktion 



X (r) = jJ d 2 r'n(r')g(r-r') 



elektrostat. Potential 



*(r) = J d 2 r> e (r')ff(r-r') 



elastische Energie 

Hel = \J d 2 rr,(r) X (r) 



Energie 



H = \J <i 2 rj(r)$(r) 



und 



V 4 X (r) = Jr,(r) 



47r»?(r) 



Poisson- Gleichung 
V 2 *(r) = - 



mit e = if. 
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Im folgenden wird die elektrodynamische Mittelung iiber atomare bzw. 
molekulare Strukturen verglichen mit der Mittelung iiber Dipole aus Burgers- 
Vektoren. Die elektrodynamische Vorgehensweise entspricht der von Jackson 
(1963). 



Elastizitatstheorie Elektrodynamik 

Die Quellenfunktio n ei ncs Dipols T ist in inneren Die Ladungsdichte des j-ten Molekiils ist in 

Koordinaten p( r ) ( fj.9| ): internen Koordinaten r': 

%V r) ) e'ii*') 

Die Verzerungsfunktion bei r erzeugt durch Das Potential bei r erzeugt durch das j-te 

Dipol T, dessen Zentrum bei r r liegt, ist Molekul mit Zentrum bei rj ist 

Xr(r) = jj d 2 p (r) *(p (r) )9(r - r r - p< r >) ^(r) = J dV ^(r') 9 (r - rj - r'). 



Entwicklung dieses Ausdrucks 



Entwicklung dieses Ausdrucks 



fur|p(r)|«| r _ rr | : ^ fur|r'|«|r- rj |: 
X r(r)= J J d 2 p^vr(p^)g(r-r T ) 

2 J J dxndxrj 



^(r) = J dV^CrOsCr-rj) 

/ ./ /■_/•> _/ d g(r-rj) 



(kein Term erster Ordnung) 

Gesamtpotential durch Summation iiber alle 
Gesamtverzerrungsfunktion durch Summation Molekiile 
iiber alle Dipole 

*(r) = ^* J -(r) 
x(r) = 2^Xr(r) j 

r 

nur 1. Term 

nur 1. Term . 

Xl(r) = j£ y d 2 p (r) »7r(p (r) )!?(r-r r ) j ^ 

/• „ = [ d 2 r'g mol (r")g(r -r") 

= J / d r 77 p (r )<?(r — r ) J 1 W ' 

mit 

mit £W(r") = Z) / dV ^'( r ')<5(r" - rj) 

»W) = £ / ^ (r) *(p (r »)i(r"-r r ) 
r 
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Mittclung iiber Flachen, die grofier sind als 
quadricrtc maximalc Dipolseparation 
AA»C 2 : 

X g(r - r") 
= jj d 2 r"( V Dip (r"))g(r-r") 



Mittelung iiber Volumina, die grofier sind als 
die maximalen Molckiilvolumina 
AV» V mol : 

<*i (>■)> = 7^7 / <i 3 e*i(r+e) 

AV J AV 

=^L d3 'J d3r " s - i{r " + ' ] 

xg(r-r") 
= J d 3 r"(e mol (r"))g(r-r") 



Damit makroskopische Qucllenfunktion Damit makroskopischc Ladungsdichtc 

(V mp (r)) = (v mp (r)) AN n(r) (e mo l(r)) = <e moi (r)> n(r) 

wobei wobei 

(f} Dip (r))^ N ■ Durchschnittsquellen- (e mo i(r)> : Durchschnittsladung 

funktion pro Dipol pro Molckiil 

n(r) : makroskop. Dipoldichte n(r) : makroskop. Molekiildichte 

d.h, wenn in AV durchschnittlich AW Molekiile: 

d.h, wenn in AA durchschnittlich AN Dipole: i r 

(e m ol(.r)) = —— d i (g mo i(r+^) 



(v Dip (r))AN = -r^ / d\r, D ^{r + i) 
AN J AA 



AN 



Erlaube nun frcic, von Molekulen unabhangigc, 
Ladungen. Die makroskopischc Ladungsdichtc 

Erlaube nun Dipole mit grofier Separation r > w i rcl zu - 

f. Die makroskopische Quellenfunktion wird ^/ r \ _ j(r)) + g ex (r) 



zu: 



7p (r) = {ri Dlp (r)) + r) e x (r) Analoge Mittelung filr zweiten Term des entwickcltcn 

PotentialsV liefert makroskopisches Dipolmoment 
Analoge Mittclung fur zweiten Term dcr entwickeftfe) • Damit wird das makroskopischc Gesamtpotcntial 
Verzerrungsfunktion liefert makroskop. Dipoltensor f 3 i r i 

C\- (r). Damit wird die makroskop. Gesamtverzerrung's^riktjorf 1 ' r )^ ( r — r ) 



x (r)=J J d 2 r' V <(r')g{r - r') + y dVpi(r')- 

+ j/d 2 r^(r0 82g(r - r,) 



9x'. 
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Partielle Integration und Ableiten liefern 

d 2 C c (i 
»j v 

dx'idx j 



V 4 x(r) c „ , /- j2 d^jr') 
__ = ^ (r) + J d r 6(r - r 



r,«(r). 



Partielle Integration und Ableiten liefern 
v 2 <£ (r) _ , f , 3 ,dpiiy). 



Air 



- C (r) + / ^'^(r-r') 



-«?( r ) + ^-P'( r ) 



Mit Linear Response geben die Symmetrien 
Mit Linear Response geben die Symmetrien deg isotropen Systems 



des Drciccksgitters 

C&(r)=£°y*(0 



a 2 x(r) 



dx k dx t 

1 c\{o( d2xir) 1 a2x(r) 

2 \dxidxj dxjdxi 



d 2 x(r) 
d 2 x k dx l 



d$(r) 



-Xei" 



9<I>(r) 

dxi 



Nun kann die dielektrische Konstante berechnet Nun kann die dielektrische Konstante berechnet 
werden werden 

£ V 2 $(r) = -47T0(r) 



£ (C)V 4 X (r) =47T^(r) 



In 



V 4 x(r) 
-47r(Ci(C) + C 2 (C))V 4 X (r) 



=V 2 *(r) + 4tt— — *(r) Xei 




Somit gilt fur e 

e = e -4tt(Ci +C* 2 ) 

Die clcktrischc Suszcptibiltat k ist also 

K=-Cl-C 2 . 



Somit gilt fur e 

£ = 1 + 4vr x e i 

Die elektrische Suszcptibiltat ist also x e i- 
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Anhang B 

Unordnungswechselwirkung 



B.l Varianz bei einer einzelnen freien Versetzung 



Entsprechend Formel ( 2.63| ) betragt die Varianz der Wechselwirkung mit der 
Unordnung unter Beriicksichtigung der Translationsinvarianz des Systems 



:(R) = [V(i 



[V(0) 



Nach (|2.54|) gilt 



V(0) = f cPrScir^j^lrl 



(B.l) 



(B.2) 



Mit Hilfe der Relationen [5c(r)Sc(r')]o = a6(r — r') und a — aVL^ ergibt sich 

'j 2 n 2 



[V(Of 



I D 



16 

J 2 5 



^bidj J d 2 rSc(r)V, In \r\b k e kl J dV<fc(r')V{ In |r'| 



D 



16tt 2 1 10 
J 2 d 



biCijbkeki I d 2 rV j In |r|Vj In |r| 



16tt 2 



hi I d 2 r V — + bl I d 2 r^-2b x b y I d 2 r^ 



(B.3) 



Unte r Ver wendung des Cut-Offs 1 fiir die Integration folgt schliefilich die Behauptung 
aus (2.6c) 



[v(oy 



J 2 a 



D 16tt 2 
J 2 



bl InR 



b 2 nlnR-0 



(B.4) 



16tt 



a InR. 
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B.2 Varianz bei einem Versetzungsdipol 



Nach (2.70) gilt fur die Varianz der Unordnungswechsclwirkung eines Versetzungdipols 

A 2 (r, 0) = [(V(r) + V(0)) 2 } D = 2 [V(0)V(0) + V(r)V(0)] D . (B.5) 

Aufgrund der Berechnungen aus B.l ist nur noch notwendig [V(r) V(0)]£> zu 
untersuchen. Mit (2.54) hat man hier bei analogem Vorgehen zu B.l 



[V(r)V(0)}, 



J 2 Q 2 



16tt s 



Xbp k e kl J d 2 r"8c{v")V'{\n\v"\ 



(B, 



D 



J 2 a 
16tt 2 



b at e tJ b( 3k e k i I d 2 r'V; in|r-r'| Vj In |r'| 



Da fiir einen Versetzungsdipol h a = — gilt, kann auf die Unterscheidung 
verzichtet werden, wenn man obige Gleichung mit (—1) multipliziert. Verwendet 
man die Relation A(r 2 In |r|) = 4 + 4 In |r|, so ergibt sich, da die Konstante bei 
nochmaliger Differentiation wegfallt 

[V{v)V{Q)) D = -\^b i e ij b k ekl I d'r'V* In |r - r'| V;A'(r' 2 In |r'|). (B.7) 



Zwcifache partielle Integration liefert als Haupterm H 

H{v,Q) = -^b i e ij b k e kl J d 2 r'V'^' In |r - r'| Vj(r' 2 In |r'| 
= -^rbfyjbkekiVjViir 2 In |rj), 

OZ7T 



(B. 



wob ei let zteres nur fiir Separationen r>l gilt. Wendet man hierauf das in ( A.l ) 
bis ( AA) beschriebene Vorgehen an, so erhalt man unter Vernachlassigung von 
konstantcn Beitragen 



i/(r,0) = -^<x(ln|r|-cos 2 

l07T 



(B.9) 



Es gilt nun die Randterme der partiellen Integration zu betrachten. In zwei 
Dimensionen hat man fiir den Laplace-Operator dargestellt in Polarkoordinatcn 
A = ^jf^rJ^: + ■pz-gfiz- Beschrankt man sich zunachst auf die erste partielle 
Integration fiir den ersten Term des Laplace-Operators, so ist der Randtermbcitrag 

R 



j2 — r o 

Rn{v,0) =-^bieijb k e k i / #r'V^ -In |r' - r| V[ — (r' 2 lnr' 



64tt 2 

Mit Hilfc der Rclationcn 



(B.10) 



(B.ll) 
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crhalt man somit 



Rri(r, 0) = -——^bieijbktki 

D47T 



J 2 



X ■ X j 

>—+ 

(r' - r) 2 



(3x[ + 2x[ In/) 



(B.12) 



32tt 



a In i?, 



wobei von x' = r'cos0 und y' = r'sin</> Gebrauch gemacht wurde. 

Die zweite partielle Integration des ersten Laplace- Terms hat den Randterm 



i? r2 (r,0) 



J 2 a 



647T 2 "* 13 
Hier gelten die Beziehungen 



-'V'^-ln|r'-r| VAr' 2 \nr') 
3 or' 



(B.13) 



V'.— lnlr'-rl w Xj X ° 

J dr' 1 1 (r'-r) 2 r' 

Vj(r' 2 lnr') = + 2a,J lnr', 
die analog zu oben eingesetzt folgendes ergeben 
J 2 a 



(B.14) 



R r2 {r,0) 



64tt 2 ^ " 
_ J*_ 
32tt 



biCijbkCki 
alnR. 



271 



' rj $ + 2^ lnr') 



(r' - r) 5 



(B.15) 



Die beiden Randterme liefern also den gleichen Beitrag. Die Beitrage durch 
die Randterme der partiellen Intergration fur den zweiten Teil des Laplace- 
Operators (p-Tjrj-) sind 0, da hier nur Integrale der Form 



drf(r, cos cj>, sin (j. 

auftreten. Es gilt also insgesamt 

[V(r)V(0)] D = H(r, 0) + R rl (r, 0) + R r2 (r, 0) 
J 2 

= -—a(lnr- cos 2 6-lnR). 

107T 



(B.16) 



(B.17) 



Abschlicficnd kann festgestellt werden, daB die Varianz fur Separationen r> 1 
den in (2.70) angegebenen Wert annimmt 



A 2 (r,9) = 2[V(0)V(0) + V(r)V(0)] D 

J2 _ 

= —a (in r — cos 2 0) . 



(B.18) 



Dieses Resultat wurde auch von Nelson g bei der Berechnung der von ihm 
benutzten replizierten Hamiltonfunktion verwendet. 
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Anhang C 



Berechnungen zur 
Polarisierbarkeit 



C.l Polarisierbarkeit in der Elektrodynamik 

Dieser Abschnitt soli zeigen, dafi die in Abschnitt 3.3 angewandte Methodik 
zur Bestimmung der Polarisierbarkeit genauso in der ublichen Elektrodynamik 
benutzt werden kann und daraus ihre Berechtigung erfahrt. 

Hierzu betrachten wir ein System von Punktladungen. Unordnung spiele keine 
Rolle. Die Wahrscheinlichkeit, an einem Paar von Platzen im System mit Abstand 
r einen Dipol zu finden, ist Boltzmann-verteilt p(r) oc exp(— 0E), E ist die 
Wechselwirkungsenergie der beiden Ladungen. 

Legt man nun ein konstantes, homogenes E-Feld an das System an, so ist die 
Wechselwirkung eines Dipols mit diesem Feld 



Hierbei ist fj, = r das Dipolmoment eines Dipols mit Separation r und Ladung 



W = -mEi. 



(C.l) 



Analog zu (3.53) laBt sich die Polarisationsdichte qi(r) als 



w + — w_ 



(C.2) 



1 + w + + w. 



mit 



w± = exp[-(3(E ± W)} 



(C.3) 



bestimmen. Auch hier kann Tensor Aij eingefiihrt werden 




(C.4) 
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Der Vorfaktor beriicksichtigt, dafi eine Orientierung in ( |C.2| ) modulo tt ausgezeichnet 
ist. Ableiten gibt 



(w + + U>_)(1 + W + + W-) — (u> + — W-)' 



E=0 



Damit kann dcr Tensor An berechnet werden 



1 Jo 2lT 

r 2 p(r) 



2T 



E=0 



(C.5) 



(C.6) 



Da in der Elektrodynamik in einem isotropen Medium die induzierte Polarisation 
parallel zum anliegenden E-Feld ist 



(C.7) 



nimmt die Polarisierbarkeitsdichte A(r), fiir die Xe( r ) — Jo dr'A(r') gilt, folgende 
Form an 



A(r) = A xx (r) = A yy (r) = 



r 2 p(r) 
2T 



(C.8) 



Die Polarisierbarkeit a(r) wird damit zu 



a(r) 



A{r) _ f_ 
p(r) ~ 2T" 



(C.9) 



dem aus dcr Elektrodynamik bekannten Ausdruck. In Abschnitt 3.3 wurde diese 
Methode zur Polarisierbarkeitsberechnung iibertragen auf das zu untersuchende 
Problem eines Coulombgases mit Vektorladungen. 



C.2 Winkelintegration des Dipoltensors 

In diesem Abschnitt soli das Integral 

^Q^Qkie*™ 26 (C.IO) 
Jo 27r 

gelost werden, welches bei der Berechnung der Polarisierbarkeitsdichte immer 
auftritt, lcdiglich der Wert von x hangt vom betrachteten Temperaturbereich 
ab. Im weiteren wird folgende Abkiirzung benutzt: 

r 2n Hf) 

<f(r,6)> w := J o — f(r,9)e xcos 9 (C.ll) 
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Da entsprechend Gleichung (3.51) gilt 



QijQkl — ^ir^j s^kt^lu{b r C srn X 7n ~t~ ^s^rm^mjiPt^un^n "f~ ^li^tn^n)) (C.12) 

beschranken wir uns zunachst auf die Berechnung von 

<C (b r c srn x m -\- b s c rm x m )(^btc un x n -)- b u €t n x n *) • (C.13) 

Mittels elementarer Integration unter Beachtung der Additionstheoreme fiir 
Cosinus und Sinus (vgl. auch Young ||l3f ) , erhalt man die folgenden Beziehung 

< XiXj > w = (^Y S v + 7,2 ^ b J - \ Si ^j cos2 ®J ■ (C14) 



Wendet man diese auf ( C.13 ) an so ergibt sich unter Verwendung von eyejw 

5ik$ij — SuSjk- 

<C (6 r c sm x m ~t- b s c rm x m )(^btc un x n -(- b u €t n x n *j ^>w 



(C.15) 



~\~ b r b u € S m€tn%m%n b s btC rm C un X m X n ^ ^>w 

— <C \brbt{S suSfan 6 sn^um) ~\~ bsbu^rt^mn $rn$tm) 
~\~ b r b u {^S s tS mr i ^sn^tm) H~ b s bt(^S ru S mn — $rn$urn)\ 

x r 2 (S mn /2 + [b m b n - S mn /2] cos26>) > w . 

Fiir die Mittelung iiber alle Richtungen des Burgers- Vektors b setzt man ein 
isotropes System an. Hier gilt 

< b m b n > m — 2^™n (C.16) 

< bkbib m b n >av= -^{Skl^mn + Skm&ln + ^kn^lm)- (C17) 

o 

Fiir das Dreiecksgitter gelten beide Gleichungen wie fur das isotrope System 
exakt, im Fall des Quadratgitters kan n (|C.17 ) als Naherung benutzt werden. 



Wendet man diese Mittelung nun auf ( |CTW an, so resultiert 



<C (pr^sm^m H~ b s € rrn X m ^){bt^un^n b u tt n X ri ^ ^>yj 

(C.18) 

2 {6 rt 5 su + 5 ru 5 ts )--co S 265 ra 8 tu 

Mit der Umformung exp(a; cos 2 9) — exp(|) exp(| cos 29) und unter Benutzung 
der modifizierten Besselfunktionen Jo und I\, die wie folgt defmiert sind 

Z7T 



(I 



(C.19) 



r27T in 

h ( x ) ■- / ™ cos2 9e xcos26 , 

Jo 27T 
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findet man schliefilich fur das gesuchte Resultat 

f 

Jo 



— n n p x cos2 9 

Wij Wkl ~ 
Z7T 



(C.20) 

1 

4' 



■je ir e js e kt ei u [Qi(r)(8 rt 5 su + 5 ru S ts ) + Q2(r)S rs 5 t u] , 



wobei 



Qi{r) = \r 2 I [^ 



I . (C21) 



Diese Relation erlaubt die Berechnung der Polarisierbarkeiten und Polarisierbarkeitsdichtcn 
in Abschnitt 3.3. 
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Anhang D 



Mittelung des 
Ordnungsparameters 

Unter Verwendung der Methoden von Nelson H soli in diesem Abschnitt der in 



Kapitel 5 ( 5.15 ) verwendete Ausdruck fur den Ordnungspar. [(exp(iGu(r)))] £) 
hergeleitet werden. Hierbei steht < ... > fur eine thermische Mittelung und 
[■■■]d fur die Mittelung iiber die moglichen Unordnungskonfigurationcn. Wir 
betrachten das System auf sehr groBen Langenskalen, wo Versetzungen in der 
kristallinen Phase kcinc Rollc mchr spielen, da sie zu Dipolen mit kleinerer 
Separation gebunden sind. Folglich kann bei der thermischen Mittelung ein rein 
phononisches Verschiebungsfeld angesetzt werden, welches Gaufi-verteilt ist. 



D.l Berechnung ohne Unordnung 

Zunachst betrachten wir ein System, das keine Verunreinigungen aufweist; in 
der Hamiltonfunktion wird Sc null gesetzt. Mit einem beliebigen reziproken 
Gittervektor G und dem Verschiebungsfeld u(r) kann man den Ordnungsparameter 
entwickcln 

(exp(*Gu(r))) . f><™ . B^M^, ( D1 ) 

n=0 ' n=0 *> >' 

die ungeraden Terme der Entwicklung fallen weg, da iiber eine gauBische Wahr- 
scheinlichkeits-Verteilung gemittelt wird. Das Wick-Theorem liefert nun 

((Gu) 2 ") = gll((Gu) 2 )». (D.2) 

Somit ist 

(exp^Gu(r))) = ^ = exp (-^((Gu) 2 )] , (D.3) 
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was es notwendig macht, den Ausdruck 

((Gu) 2 ) - Gl(ul) + 2G x G y (u x u y ) + G 2 y (u 2 y ) (D.4) 

zu berechnen. 

Urn die Hamiltonfmiktion 

H=\J d 2 r2fiu 2 k (r) + \u% (r) (D.5) 

im Impulsraum darstellen zu konnen, fiihren wir einen Basiswechsel von einer 
Darstellung (u x ,u y ) in eine Darstellung (ui,u t ) durch. Hierbei ist ui(k) der 
longitudinale Anteil von u beziiglich k und Ut (k) der transversale Anteil 

ui = e k (e k u) (D.6) 
u t = u ui = u e k (e k u). (D.7) 

e k symbolisiert den Einheitsvektor in k-Richtung. Die zugehorige Basiswechselmatrix 
ist 

A=y^ f i, (d. 

und die zugehorige inverse Matrix 



fcx 


ky 


VP 


VP 


ky 


k x 


VP 


VP , 


k x 


ky 


VP 


~ VP 




fern 


VP 


VP 



(D.9) 

Mit diesem Basiswechsel wird die fouriertransformierte Hamiltonfmiktion zu 
(vgl. auch 

H=\J ^{2^ + \)k 2 u 2 + ^k 2 u 2 (D.10) 



Durch Ableiten hndet man 

dH 



= {X + 2n)k 2 ui(k) (D.ll) 



du t {k) 

Benutzung des Aquipartitionsprinzips 

dH 

liefert 



r) H 

Hk 2 u t (k). (D.12) 



u i— )= T§ ij (D.13) 



und 



1 c) M T 

= (uW <Ul W = (XTW (m4) 



(u t (k)u t (k)) = -j^. (D.15) 
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Ftihrt man jetzt die Riicktransformation mit Hilfc der inversen Matrix A 1 aus, 
so resultiert unter Beachtung von (u t (k)u/(k)) = (wegen (13.13)): 



(Ui(k)uj(k)) = fc l kj{uiui) + (Sij - h k.j)(u t u t ) 



T 



T 

fa kj + ~~^2 ~ kj ) . 



(D.16) 



Hierbei stehen die fur Kombinationen der (x,y), ferner ist fcj die z-te 

Komponente des Einheitsvektors e^. Fiihrt man nun eine Fouriertransformation 
durch, so erhalt man wegen J Q n ki kjdO = n5ij folgenden Ausdruck 



i( T )uj{r)) = 



1 



(2tt)2 
13 Air 



d 2 fc<Ui(k)uj(k)) 
T 



dk k 



1 



= T hxR 



hxR 



(A + 2n)k 2 
T T 



T 
Jik 2 



(A + 2[i fj, 
3fi + X 



47T/i(A + 2ju) 
(K(r)) 2 ) = (K(r)) 2 ). 



Setzt man dieses Resultat in ( DA ) ein, so findet man 



((Gu) 2 ) = G 2 T\tlR 



3// + A 



47r/x(A + ' 



(D.17) 



(D.18) 



Mit (D.3) erhalt man schlieBlich im Fall ohne Unordnung das bereits von Pokrovsky 
p6[ gefundene Ergebnis 



T G (3j£+X) 

(exp(iGu(r))) = i? 8 -f( A + 2 ^) . 



(D.19) 



D.2 Berechnung mit Unordnung 



Unter Beriicksichtigung der Unordnung nimmt die zu betrachtende Hamiltonfunktion 
die folgende Gestalt an 



H = H - J d 2 r(fi + X)fl Sc(v)u kk (r). 



(D.20) 



Dabei stent Hq fiir die Hamiltonfunktion des Systems ohne Unordnung (D.5). 
Ebenso bezeichnen < ... > und < ... >o im folgenden die thermische Mittelung 
iiber den vollen Hamiltonian, bczichungsweise die Mittelung mit Hq. Damit ist 



[(exp(iGu(r))}] fl = 
Setzt man nun 



J G u(r) e (A1+ y n ° / d 2 rfcVuj 



D 



T 



d rScVu, 



(D.21) 



(D.22) 
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so gilt 

(e*)„ = e M x ^°, (d.23) 
da die Mittelung iiber Hq wiederum eine gaufiische Mittelung ist. Wir finden 



} =-<(Gu) 2 ) + 2 



T 



. Q + A)O 
1 T 



d 2 rScGi ( Ui (0) V j u J -(r)) 



(D.24) 



d 2 r5c\J\i)l. 



Setzt man dies in (D.21) ein, so heben sich Nenner und der letzte Term des 
obigen Ausdrucks weg. Es bleibt: 



[<e X p(iGu(r)))] B 



-i((Gu) 2 ) + g » U ' + y n ° J d 2 rScG i V j {u i (Q)u j (r)) 



D 



-i((Gu) 2 ) +e t/(0) 



(D.25) 



Die Unordmmgsmittelung des ersten Terms ist unnStig, da er nicht von der 
Unordnung abhangt, er wurde in Abschnitt D.l berechnet, das Resultat ist 



( D . 1 9 ) . Da die Unordnung gauBisch verteilt ist , gilt fur die Unordnungsmittelung 
des zweiten Terms 



J 1(0) 



D 



(D.26) 



mit [Sc(r)Sc(r')] D — aS(r — r'). Daraus folgt 



(m + A)n 

T 



aG r G s / <2 2 rVi(u r (G) tti (r)) V^u^O) w i (r)) . 



(D.27) 

Mit Hilfe der in D.l beschriebenen Methodik erhalt man hier folgende Darstellung 



(2/x + A) 2 



G r G s 



(2tt) 2 fc 4 



(M + A) 2 » 2 a 
47r(2/i + A) 2 



(G 2 +G 2 ) / dkk 



1 

'fc 2 



(D.28) 



lni? 



G 2 a(fi + X) 2 n 2 a 
4tt(2^ + A) 2 ' 



Setzt man dies in ( |D.25| ) ein, so ergibt sich die in Kapitel 5 benutzte Formel 

(HI) 



[(exp(iGu(r)))] D = R *^m^l R (D.29) 
die in ahnlicher Form bereits von Nelson [^) hergeleitet worden ist. 
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